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ABSTRACT 


The Single Stage Capacitated Warehouse Location Problem (SSCWLP) has been 
formulated in two different Styles as given by Geoffrion and Graves (1974) and Sharma 
(1991). We have used the two styles and combined them with the ideas of strong and 
weak relaxations ( as picked up from the formulations of capacitated plant location 
problems) to develop variety of new formulations of SSCWLP. Later we also added 
Big-M constraints that link the real variables (distribution variables) with the 0-1 integer 
variables (location variables). 

We empirically established the strengths of Linear Programming relaxation of various 
formulations developed above. Among significant findings we found that Big-M 
formulations significantly boost the bounds (over and above the bounds given by WRS 
and SRS) for the case when Ecapj= 2.0 and the constraint cap J y J >=Z,i c xp. However it is 
not true when Scapj= 6.0 or capj>=E, k Xjj k 

Another Significant finding is that for solving linear programming relaxations of all 
kinds, the style of Geoffrion and Graves (1974) took significantly more computational 
iterations than the style of Sharma (1991). 

We also found that the constraint cap J y J >=E, k x;j k improved the bounds substantially in all 


possible formulations. 
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Chapter 1 


INTRODUCTION 


A commonly occurring problem in distribution system design is the optimal location of 
intermediate distribution facilities or the warehouses between plants and customers. A 
capacitated single period version of this problem is formulated as a mixed integer linear 
programming problem in differing styles in literature. 

The formal definition of the problem considered in this work and referred as a Single 
Stage Capacitated Warehouse location problem can be given as below. The problem 
SSCWLP arises when the distances between plants and markets are large and it 
becomes necessary to route the supplies through warehouses known as trans-shipment 
points. The set of potential trans-shipments locations are known and each point has an 
associated fixed cost with it The problem is to chose a sufficient number of trans- 
shipment points such that the sum total of fixed location costs and transportation cost of 
shipping goods to market is (from plant to warehouse and from warehouse to market) 
minimized while satisfying the demands at each market points. It is assumed that the 
warehouses have finite capacity and which is known quantity. Also only a single 
commodity is considered for distribution. 

Or in other general facility location is the problem of locating a number of facilities 
from a subset of m potential facility locations to completely satisfy at minimum cost all 
the demands of a set of n clients. Each client has an associated demand and the costs 
considered include transportation and fixed costs for opening facilities. For the 
Uncapacitated Warehouse Location Problem ( UWLP ) the facilities are assumed to have 



unlimited capacity, and for the Capacitated Warehouse Location Problem ( CWLP ) 
there are constraints on the total demand that can be met from a facility. 

Both the problems are well known to be NP-hard. 

That is, if the number of candidate warehouses (m) and customers (n) are growing 
larger, the number of constraint equations for Integer Programming explodes and it is 
very hard not only to solve the problem, but even to represent it in the memory of a 
computer. Therefore, we need a practical method for finding an approximately optimal 
number and locations from a large number of possible locations. 

So a multicommodity warehouse location problem may have a lot of aspects to consider 
such as supplying costs, supplying limits, list of commodities, list of eligible warehouse 
locations, minimum and maximum limits of throughput for the facilities, fixed cost of 
setting up the warehouse, customer demands, transportation or freight costs etc. 

These considerations add up to the complexity of the model so a practical method is 
needed for real life large sized problems so that sufficiently closer bounds can be 
obtained with relatively less computational efforts. 

To explain the requirement and function of the solver the diagrammatic sketch of a 
warehouse location problem is given as below. 



FIGURE 1. SKETCH OF A COMPREHENSIVE DISTRIBUTION 
PLANNING SYSTEM. 
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MAIN FUNCTION OF THE SOLVER 


In this work the Single Stage Capacitated Warehouse Location Problem (SSCWLP) is 
formulated by two different styles. In style 1 the commodity is considered to be 
distributed from the supply locations to the markets in one stage. Style 1 is given in 
literature by Geoffrion and Graves (1974) [1], In style 2 two stages are considered once 
the commodity goes from the plants to the warehouse and in the second stage from the 
warehouse to the markets. This style is given in literature by Sharma (1991) [2]. For 
each of the two given styles many relaxations are given and is emperically seen that a 
new relaxation (Big M method) gives a better bound than all the relaxations tried. 
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The organization of the thesis work is as follows: 


In chapter two, we have given the relevant literature review. In chapter three we list 
down all the formulations for the SSCWLP for both the styles. In chapter four the 
computational experience and results for problems of various sizes are solved. In 
chapter five the analysis of the data in chapter four is done to establish which relaxation 
performs better and when and also the conclusion is given. In chapter six we have 
included a new formulation after getting good results from one of our formulations and 
the computational experience and analysis of the data for this formulation. Chapter 
seven gives the direction for future research on the topic of interest. 

The formulations have been solved by using the optimization package LINGO and the 
objective values, and no. of iterations needed are noted for the various formulations. 
Approximately 2200 problems of various sizes were solved and each problem was 
randomly generated for this purpose. 
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Chapter 2 


LITERATURE REVIEW 

In this chapter, we give a brief literature review on SSCWLP (Single Stage Capacitated 
Warehouse Location Problem). 

2.1 Early Heuristics 

One of the earliest methods proposed for warehouse location/ facility location is the 
now well known heuristics proposed by Kuhen & Hamburger (1962) [13]. It has been 
widely studied and its catalytic effect can hardly be overrated. The heuristic consists of 
two parts. The main program is an ‘add routine’ whereby facilities are located one by 
one corresponding to the greatest cost reduction until no facility can be added without 
increasing the total cost. The underlying hypothesis is that the optimal solution for 
‘(p+1)’ facilities can be determined from the optimal solution for ‘p’ facilities by adding 
an additional facility to the existing solution. In the modem technology, such a scheme 
would be called greedy algorithm because of its appetite for maximum improvement at 
each step. Upon termination of the main program, a so called bump and shift routine is 
entered. It first eliminates (bumps) any facility which is now uneconomical because of 
the proximity of another facility located subsequently. It also considers relocating 
(shifting) each facility from its actual location to other potential locations in its 
neighborhood. This was first applied to simple plant location problem; and possibly can 
be applied to SSCWLP. 
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2.2 Genetic Algorithm 

Genetic Algorithms (GA’s) are robust and adaptive methods that may be used to solve 
the problems related to search and optimization. GA’s work with a population of 
individuals (usually 10-200), each representing a possible solution for the given 
problem. To each individual fitness value has assigned according to adaptation of this 
individual. The population evolves towards better solutions by means of randomized 
processes of selection, crossover, and mutation. The selection mechanism favors 
individuals of better fitness value to reproduce more often than the worse ones when 
forming a new population. The crossover allows mixing of parental information when it 
is passed to their descendants. The result of crossover is structured and randomized 
exchange of genetic material between solutions, with the possibility that "good" 
solutions can generate "better" ones. The mutation changes a bit in the binary string 
with some small probability pm. The role of mutation in GA’s is restoring lost or 
unexplored genetic material into the population. It can be used to prevent the premature 
convergence of GA to suboptimal solutions. The initial population is usually randomly 
initialized. 

Kratica, Filipovice and Tosice (1998) [6] have used a new hybrid approach for solving 
the Warehouse Location Problem by Simple Genetic Algorithm (SGA) and Add- 
Heuristic. The best individual in every generation of SGA is improved by Add- 
Heuristic, adding or removing one warehouse with maximal reduction of overall cost. 
The existence of warehouse with adding/removing may reduce overall cost and binary 
string is changed in the corresponding position. SGA itself produces good results with 
the reasonable running time, but hybrid approach improves performance of 
implementation. The local search nature of Add-Heuristic provides achieving better 
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regions of search for given problem. The GA implementation assisted with Add- 
Heuristic successfully prevents possibility of premature convergence, and restores lost 
or unexplored genetic material into the population. During the testing process, in both 
cases (base GA and hybrid GA) for every particular test problem at least one-half of 
optimal solutions are obtained. Both implementations are able to solve large-scale 
warehouse location problems, with good running time. According to WLP benchmark, 
the use of hybrid GA technique brings some improvements. In general, the average run- 
time produced by hybrid GA implementation is smaller than the one produced by the 
base GA implementation. The quality of solutions and average number of generations 
are also better for hybrid GA implementation. 

For the given problem binary coding of possible warehouse locations is used. Every bit 
in the binary string denotes that particular warehouse is established if its value is 1, 
while 0 denotes it is not established. For the faster execution the item string is allocated 
in 32-bit words. The objective value function is computed depends on number of 
established warehouses e and threshold eO. If e > eO in initialization part, warehouses 
are indexed in non decreasing order of transportation costs for every customer. After 
that the objective value function speedly finds first established warehouse, with minimal 
transportation cost. If e <= eO previous procedure is not appropriate, and ordinal array 
of established warehouses is formed. We find the established warehouse with minimal 
cost searching only the same ordinal array for every customer. The whole population is 
replaced in every generation except the best individual. The best individual in the 
current generation directly goes to the next generation, without selection, crossover and 
mutation. The initial population is randomly initialized, because the maximal diversity 
into the population is maintained in that case. The duplicate strings are discarded from 
the population, i.e. multiple occurrence of the same string is practically excluded. The 
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duplicate strings in every generation are discarded by setting their fitnesse’s to zero. 
This method does not remove duplicate item string physically, but eliminates its 
occurrence in the next generation. This technique maintains the diversity of genetic 
material. Discarding of duplicate strings decreases the appearance of dominating 
individuals in the population, and effectively decreases the possibility of premature 
convergence in local optimum. This is a very important factor for successful working 
GA, particularly in the cases of large size test problems. 


23 Benders Decomposition 

Geoffrion and Graves (1974) [1] one of the earliest researchers in the area of interest 
have formulated a multicommodity capacitated version of the problem as a mixed 
integer linear program. A solution procedure based on Bender’s Decomposition is 
developed, implemented and successfully applied to a real problem for a major food 
firm. An essentially optimal solution was found and proven with a surprisingly small 
number of Bender’s cuts. The major conclusions arising from the study is the 
remarkable effectiveness of Bender’s Decomposition as a computational strategy for 
static multicommodity intermediate location problems. The numerical experience 
quoted also shows that only a few cuts are needed to find and verify the solution within 
one or two tenths of the global optimum. 

The mathematical formulation of the problem uses the following notation 

T Index for commodities 

‘j 5 Index for plants 

c k’ Index for possible distribution centre or warehouse 
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‘ 1 ’ 

Sij 


Du 


Vk,Vk 


fi 


Vk 

Cijkt 


Xijkl 


Index for customer demand zones 

Supply (production capacity) for commodity T at plant ‘j’ 

Demand for commodity T in customer zone T 

Maximum and minimum allowed total annual throughput for warehouse 
at site ‘k’ 

Fixed set up and operating cost for the warehouse at site ’k’ 

Variable unit cost of throughput for a warehouse at site ‘k’ 

Average unit cost of producing and shipping commodity ‘i’ from plant 
‘j’ through warehouse ‘k’ to customer zone T 
A variable denoting the amount of commodity T shipped from plant ‘j 5 
through warehouse ; k’ to customer zone ‘1’ 


A 0-1 variable that will be 1 if warehouse ‘k’ serves customer zone T 0 
otherwise 

A 0-1 variable that will be 1 if a warehouse is acquired at ‘k’ 0 otherwise 


) Xijkl.Cijkl + y 

MINIMIZE jjti k 

SUBJECT TO- 


^ * Xijkl — SijXf jj ( 1 ) 

k,l 


fkZk 


+ V/t^T Duyki 

n 


^ xijkl = D uykN ikl 


j 



k 


(3) 


( 2 ) 
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(4) 


VkZk < 'y'Dityid < VkZkik 


il 


Linear configuration constraints on y and/or z 


( 5 ) 


Constraint 1 being the supply constraint, constraint 2 stipulates both that the legitimate 
demand must be met (when yki= 1) and that Xijki must be 0 for all ij when yki= 
0. Constraint 3 specifies that each customer zone must be served by a single warehouse. 
Constraint 4 enforces the correct logical relationship between y and z and also keeps the 
total annual throughput between given maximum and minimum limits. 

The author says that instead of using two sets of triply subscripted variables in two 
stages i.e. one from plant location to warehouse and other from warehouse to customer 
zones in the next stage, to use a quadruple subscripted formulation gives more 
flexibility. As the former suffers from the lack of flexibility for some applications 
because it “forgets” the origin of a commodity once it arrives at the warehouse which 
can be a serious drawback in real applications. Another advantage of the quadrupled 
subscripted formulation over the other style arises when some commodities are 
perishable. In such cases it may be necessary to disallow the possibility of shipping such 
commodities over jkl routes for which the total journey times are excessive. It also 
makes easy to accommodate direct plant-customer zone shipments in the proposed 
method. 

Another unique feature of the proposed model is that no customer zone is allowed to 
deal with more than one warehouse, since yj/s must be 0 or 1 and not fractional. Thus 
each customers’ demand must be satisfied by a single warehouse or directly from the 
plant. 

Since most real life applications for the proposed formulation are too large at that to be 
solved economically by existing mixed integer programming codes the authors have 
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suggested the technique of Bender’s Decomposition as the technique offers the 
possibility of making sequences of related runs in considerably reduced computing 
times as compared doing each run independently. The most striking conclusion from the 
wok is that it can be seen from the results of the problems solved by using the proposed 
model. It is seen that small number of iterations are required for convergence even with 
very small values of optimality tolerances. 

Sharma (1991) [2] has also used Benders’ decomposition as a solution approach to 
solve a fertilizer production - distribution problem in Indian context; however it was 
used along with Lagrangean method also. According to the scope of the research work 
in the real life circumstances although the rail network in India is quite large, is not able 
to directly reach most of the fertilizer consumption centers. Therefore, in a typical 
distribution system the fertilizers are moved fro the manufacturing plants to certain 
points known as rake points. From these rake points fertilizers are moved by road using 
heavy commercial vehicles to secondary points and finally reach the markets in light 
commercial vehicles. So essentially the formulation proposed is a two stage 
warehousing formulation with inventories also one of the factors included in the cost 
function. The decision making have to include in multi period context, at which rake 
and secondary points to locate warehouses, the size of the warehouses, the amount of 
inventory to be stocked at these warehouses, the transportation plan to meet the demand 
and how much to produce at each plant in each sub period. Severe restrictions exist on 
the quantities that can be transported in various directions as the number of rail wagons 
and HCV’s are limited. The warehouse space available at various points is also limited. 
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In this paper the plants, rake points, secondaiy points and markets, will be referred to as 
points of type ‘ 1 ’,’2’, ’3’, ’4’ respectively. 

The mathematical formulation of the above problem given by the author is give as 
below: 

Constants of the Problem 

The cropping season of six months is divided into K number of periods (assumed even); 
H is the number of products and U denotes the number of nutrients. The quantity of 
nutrient demanded of type ‘u’ at market ‘j’( point of type 4) in period ‘k’ is denoted by 
D^icu while fhu is the fraction of nutrient of type ‘u’ obtainable from product of type ‘h’. 
The total number of points of category T is N (i); SN (i, j) is the set of points of type 
i+1 to which a point of type ‘i’ having a number ‘j’ can supply a fertilizer, RN (i, j) 
denotes the set of points of type i-1 from which a point of type T having a number ‘j’ 
can receive fertilizer. Variables defined above have lower limit, upper limit and 
associated cost- coefficient and these constants are denoted by putting ‘L% ‘U’ and ‘C’ 
respectively before their names. 

Distribution cost is composed of production cost: CPijkh* Pijkh; cost of space 
CSljj*Slij+CS2jj*S2jj; cost of carrying inventory 0.5 * CIijkh*(Eijkh+Bijkh); and the cost 
of locating a warehouse is CLy*Lij. 

Variable Defijiition 

The beginning and the ending inventory at point number ‘j’ of point type T in the 
period ‘k’ of product type ‘h’ is represented by Bjjkh, Eyi* ; Ly is the location variable 
which is 1 if it is decided to locate a warehouse at point number "j’ of point type ‘i’ (i= 
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2 to 3) and 0 otherwise. The quantity produced at plant number ‘j’(point of type 1) in 
period k of product type ‘h’ is denoted by Puja,; QTR^h, is the quantity received at 
market ‘j’(point of type 4) in period ‘k’ of product type ‘h’. The space booked at point 
number ‘j’ of point type ‘i’ in the first and last three months of a six monthly season is 
denoted by Sly, S2ij respectively. The total quantity transported from point number jl 
of point type il to point number j2 of point type i2 in period ‘k’ of product type ‘h’ is 
represented by T.iji ;7 j7n, : TTyic denotes the total quantity transported from point 
number ‘j’ of point type T in period ‘k\ 

So the objective function takes the following form including all the above components: 


MINIMIZE 

) CP ljkh* P ljkh + * CTijkpmn * Tijkpmn + ^ ' 0.5 * CIijlch*(EijkhT'Bijkh) 

jkh ijkpmn ijkh 

+£ CSlij*Slij+^ CS2ij*S2i j + £CLij*Lij 

U ij V 

SUBJECT TO 

Subproblem constraints 

P 1 jkh + B\ jkh — E\ jkh + ^ ^ T 1 jlmkh 

m<=SN(ljj (6) 

for all j= 1 , ,N(1); k= 1, ? K ; h= 1, H 

T (i - I )mijkh + Bijkh — Eijkh + Tij(i + X)mkh 

meRN(iJ) meSN(iJ) (7) 

for all T= j= 1, N(i) ; k= 1, K and h= 1, H 

Lower and Upper Limit Constraints on Variables Defined (8) 



13 



QTRa jkh — ^ T 3mA jkh ( 9 ) 

meRN(AJ) 

For all j,k,h 

For the propose model there are K*H number of sub problems and each sub problems 

have the constraints 1 to 4. 1 and 2 are material balance conditions at plant, rake points 

and the secondary storage points. 3 is to set lower and upper bounds on variables 

whereas 4 sees that demand is satisfied. 

Linking Constraints 
H 


E fi- * Q TR a jkh > Da jku 

h = i 

For all u= 1,....U; j= 1, M; k= 1, K 


( 10 ) 


Bij(k + 1 )h — Eijkh 


( 11 ) 


For i= 1,2, 3; j= 1, ,N(i); k= 1, K-landh=l, H 


H 

TTijk = E E Tij(i + 1 )xkh 

h=l x<=SN(i,j) 

For i= 1,2,3; j= 1, N(i) ; k= 1, K 

H 

Eijkh <Shj 

h = 1 

Fori= 1,2,3; j= 1, N(i); k= l,....K/2 

H 

Y Eijkh < S 2 ij ( 14 ) 

h = 1 

For i= 1,2,3; j= 1, N(i); k=K/2+l, K 


( 12 ) 


(13) 


Location constraint 
Ur (°> 1) for all i, j 


( 15 ) 
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Constraint (16) gives a set of linear constraints representing the condition that if the 
warehouse is located at a particular point then the space booked must be within the 
permissible upper and lower limits 

Constraints (17) are the non negativity constraints of all real variables 

The author has used the application of Bender’s Decomposition to solve the above 

defined problem. The primal problem was attempted by two Lagrangean relaxation 

procedures. 

The use of the two relaxations is given as below in which a single product and a single 
nutrient is considered for purpose of illustration. 


Relaxation 1:LR1 

The constraints of the type Bx= b in the above formulation are added to the objective 
function with the multipliers. 


Bij(k + 1 )h — Eijkh 

For i= 1, 2, 3; j= 1, , N(i); k= 1, K-l and h= 1, H 


(ID 


H 


^ * Eijkh ^ iSly 


h - 1 


Fori= 1,2,3; j= 1, N(i); k= l,....K/2 


H 

^ * Eijkh 5; S2ij 


/i=l 


(13) 


(14) 


For i= 1,2,3; j= 1, N(i); k=K/2+l, K 

Multipliers for the equations 6, 8, 9 are uyk, vyk, wy k respectively. Starting values of 

these multipliers are chosen as follows 

Ufjk= 0 for i= 1, 2, 3 ; j= 1, N(i); k= 1, K-l 

Vijk, Wijk = min (CSljj, CS2y) for i= 2,3 ; j= 1, N(i); and all k 
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Relaxation 2:LR2 

New constraints were prepared as follows. 


(ii) 

H 


(19) 


Multipliers for the above equations are denoted by up, vlp, v2p, wlp, w2p 
respectively. 

Starting values for the multipliers were chosen as follows: 
up= 0 for i= 1,2, 3 ; j= 1, N(i); k= 1, K-l 

vlp, wl jj ki v2i jk , w2jj k = 0.5 * min (CSly, CS2jj) for i= 2,3 ; j= 1, N(i); and all k 

The Lagrangean relaxation procedure is stopped when the sum of total booked 
warehouse space attains the minimum level to be able to meet the demand. A backward 
recursive procedure is then used to get a feasible primal solution. That is, for the relative 
cost coefficients that are passed down to the sub problems, the last period (K ) is solved 
first and its beginning inventories are passed down to previous period (K-l) problem as 


Rij(k + l)h — Eijkh 

For i= 1,2, 3; j= 1, 


H 

N, Eijkh < S\ij 

/;= 1 

For i= 1,2,3; j= 1, 

H 

Eijkh < S 2 ij 

h = 1 

For i= 1,2,3; j= 1, 
B ij (k + 

For i= 1,2,3; j= 1, 

Bij(k + 

For i= 1,2,3; j= 1, 


,N(i); k=l. 


,N(i); k=l,....K/2 


N(i); k= K/2+1, 


N(i); k=K/2, 


K-l and h= 1,.. 

(13) 

(14) 

.K 

(18) 
KJ2 

K-l 



1 )/ < S If/' 

■N(i); k= 1, 

1 ) / ^ S 2 ij 
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ending inventories, and then problem of period (K-l) is solved, and again its beginning 
inventories are passed down to period (K-2). This recursive procedure is followed till 
the first problem is solved. Finally, the warehouse space booked is computed by setting 
it to the maximum of the ending inventory. 

It was found that, on an average for problems of different sizes LR2 converged in 
around one third the number of iterations taken by LR1. In LR2 the relative cost 
coefficient of warehouse space booked is a function of multipliers which determine the 
relative cost coefficient of beginning and ending inventories, whereas in LR1 the 
relative cost coefficient of the warehouse space booked is a function of multipliers that 
determine the relative cost coefficient of the ending inventories only. This could result 
in the better performance of LR2 than LR1. 

So the model developed here serves as a good aid to making decisions and routine 
sensitivity analysis is possible as complete information on dual variables is available. 

2.4 Lagrangean Relaxation 

Sharma (1996) [4] presents a mixed integer linear programming formulation as faced by 
the Food Corporation of India (FCI) which turns out to be single stage warehouse 
location problem (SSWLP), and suggest a branch and bound based procedure for the 
solution. At each node the branch and bound procedure results a minimum cost flow 
problem which can be solved in O (n 2 In (n)) time by using best known procedures. The 
paper presents a Lagrangean based heuristic procedure which obtains a good bound for 
the associated minimum cost flow problem in 0(n 2 ) time thus offering substantial 
saving in computations which become particularly significant while solving a large 
sized real life problem as faced by the FCI. 
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After giving a mixed 0-1 integer linear programming formulation of the single stage 
warehouse location problem the author relaxes the integer constraints in the problem to 
obtain the linear programming relaxation of the problem. In the linear programming 
relaxation of the problem further few constraints are relaxed and it is shown that the 
reduced problem can be solved optimally by an efficient heuristic. Also profitable use of 
the same method is used in Lagrangean based solution procedure to obtain a good 
bound to the optimal solution value of the linear programming relaxation. The optimal 
integer solution to the problem is obtained by initiating a branch and bound solution 
procedure. 

Problem Formulation 

In the formulation plants, trans-shipment points and the markets are denoted by points 
of type 5 iy 2’ and ‘3 5 respectively 

Constants of the problem 

Sj is the maximum quantity that can be supplied by a plant ‘j 5 ; Dk is the demand at 
market 4 k 5 ; fj is the fixed cost of locating a trans-shipment point; cy is the cost of 
transporting TDk units of goods from point number c j 5 of type 1 to point number T of 
type 2; Bik is the cost of transporting ZDk units of goods from point number T of type 2 
to point number c k 5 of type 3. 

Variable Definition 

Yy is the quantity received (as a fraction of the total demand) by the point number T of 
type 2 from point number 5 j 5 of type 1 ; Y,* is the location variable which is 1 if it is 
decided to locate a warehouse at point number T of type 2, 0 otherwise. Xjk is the 


18 



fraction of demand at point number ‘k’ of type 3 which is supplied by point ‘i’ of type 

2. 


MINIMIZE 


V" 1 " C 1 ’T * "! "C h s — * * * 


SUBJECT TO 

E2>=i 

j i 

£ Yu < Sj /(£ D*)Vj 

i k 

Yi > ]T Yi X i 

j 

]T Xtk = IV k 

i 

Xik < 7/V z, k 


( 20 ) 

( 21 ) 

( 22 ) 

(23) 

(24) 


2] Yij = 2] Xik(Dk / ^ Dk)\/i (25) 

7 & £ 

Yij>0Vi 9 j (26) 

Xik > 0 V z, k ( 27 ) 

7/ = (0,1) V Z (28) 


Constraint 20 ensures that adequate quantities are shipped from plants to warehouses so 
that total demand can be met. Equation 21 allows plants to supply within their 
capacities, 22 ensures that a warehouse is located only when the quantity shipped to the 
warehouse point is greater than zero. Equation 23 is for meeting the demand at each 
market point and 24 ensures that the warehouse is located only when the outflow from 
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the same is positive. Equation 25 are the flow balance constraints at each warehouse 
point i.e. inflow at a warehouse point from the plants is equal to the outflow to the 
markets.26 and 27 are non negativity constraint whereas 28 is the integer constraint. 

Next the author relaxes constraint 25 and the resulting linear programming formulation 
is as follows. 


MINIMIZE 

j i i k i 

Subject to 

20 to 28 and Y;>= 0 for all ‘i’. 

The new problem is clearly is the minimum cost flow network problem for which best 
known polynomial algorithms 0(n 2 ln(n)) are available. Now the author associates 
multipliers with equations 24 and 25 and include them in the objective function. And 
then he goes on to show that the reduced problem is solvable in O (n 2 ) time. 

The main problem so gets decomposed into two problems as below: 


MINIMIZE 


£ 




j L / 


+/*£(/>- 2>)k 


Subject to 
20 to 22 and 26 


Where 0<= f<= 1 , f is a constant and known 


(29) 


And 
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MINIMIZE 


(I-/)*E(/i-I>> K+ ( 30 ) 

i k 

^ 1 Xik(^B ik + ?uik — JUikJDk / ^ ^ Dk 

i k 

Subject to 
23 and 27 

Where 0<= f<= 1, f is a constant and known 

Fraction f of the fixed costs is included in the first problem whereas the remaining 
portion of the fixed costs is included in the second part. 

The author then gives heuristics to solve the two parts of the problem and it is shown 
that the optimal solutions for these problems can be obtained in O (n 2 ) time. 

2.5 Simulation Based Approach 

Hidaka , Okano (1997) [5] have used simulation for large-scale Warehouse (Facility) 
Location Problems in the real world on a digital map, and found an approximate 
solution. The problem was to find a near-optimal solution for the number and locations 
of warehouses that minimize the sum of the transportation cost and fixed cost, and meet 
the needs of 6,800 customers. The network data of the digital map were efficiently used 
to obtain candidate warehouse locations, to simulate the transportation cost, to simulate 
the warehouse fixed cost, and to find a near-optimal solution for the number and 
locations of warehouses. 

In order to solve a real instance of a large-scale warehouse location, they developed 
simulation-based methods. That is, they simulated the warehouse location problem in 
the real world on a digital map. 
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First, all warehouses and customers were assumed to be located at nodes of a network 
covering the whole given area. From the customer nodes, sets of candidate warehouse 
locations were selected, and the fixed and transportation costs for each candidate were 
calculated on the basis of the real data for the current 1 1 warehouses. The transportation 
cost was assumed to be proportional to the delivery time. 

Next, they optimized the number and locations of the warehouses in two steps. The first 
step was to apply Greedy - Interchange heuristics with the warehouse candidates 
selected in the above manner. Using these heuristics, they found an approximately 
optimal solution for the number and locations of warehouses selected from the 
warehouse candidates. 

Next, in order to improve the solution obtained in the first step, they devised a heuristic 
procedure for finding medians, and named it “Balloon Search." This procedure may be 
outlined as follows: they relocated each warehouse obtained in the first-step solution to 
an adjacent node of the network that was not selected as a warehouse candidate node. 
Warehouse T was relocated by aiming to find one median of the corresponding sub 
network that included only the subset of customers who were assigned to warehouse ’i\ 
In this process, they first reduced the subset of customers so that it included only 
customers whose transportation costs were relatively low. Then, they increased the 
number of the customers in the subset step by step, and repeatedly found medians. 
During these processes, they assumed that the fixed cost of the adjacent node was the 
same as that of the warehouse obtained in the first-step. After the first and second 
optimization steps, a near-optimal solution for the number and locations of warehouses 
was expected to be found from the large number of network nodes of the digital map. 
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2.6 Discussion 


Geoffrion and Graves (1974) [1] and Sharma (1991) [2] have given different styles of 
formulations of SSCWLP. We use these formulation styles and add a variety of linking 
constraints (that link binary 0-1 variable or location variable to real variables that denote 
quantity shipped from plant to market through warehouse). This results in variety of 
different formulations of SSCWLP. We relax the integer 0-1 constraint in these 
formulations to obtain various relaxations. These are given in chapter 3. Later we 
carryout an empirical investigation to determine strengths of various relaxations. This is 
done in chapter 4. 
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Chapter 3 


New formulations for the SSCWLP 


In this chapter we will enlist all the new formulations of SSCWLP. In the Chapter 2 
during the literature review we see that there are two major styles of mathematical 
model formulation for the SSWLP. In the first style the whole distribution of 
commodity from plant to markets is covered in a single phase as done by Geoffrion and 
Graves (1974) [1], The other style is a variation of this as used by Sharma (1991) [2] in 
which he uses a two phase for distribution of commodities (first from the plant to the 
warehouse and then from the warehouse to the markets). 

50 in this chapter also we have used both the styles and called the first style as STYLE 1 
and the latter as STYLE2.In the following section we will list our formulations for both 
the styles one by one with all the relaxations. In each style there are two versions of all 
the formulations which are given in the chapter. For the purpose of illustration; only one 
commodity is considered. 

3.1 Problem formulation (STYLE 1 and Variation 1) 

In this section we will enlist the various formulations of style 1 as adopted by Geoffrion 
and Graves (1974) [1], 

3.1.1 Constants Definition 

D k Demand for the commodity at market ‘k’ 

d k D k /SDk demand at market ‘k’ as a fraction of total market demand. 

51 Supply available at plant ‘i’ of the commodity 

Si Sj/EDk supply available at plant ‘i’ as a fraction of the total market 

demand 
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f] Fixed cost of locating a warehouse at ‘j ’ 

Cijk Cost of transporting LDk quantity of goods from T to ‘j’ to market ‘k’ 

CAPj Capacity of a warehouse ‘j’ 

capj CAPj/EDk capacity of a warehouse at ‘j’ as a fraction of the total market 

demand 

3.1.2 Variable Definition 

Xyk quantity of commodity transported from plant ‘i’, to warehouse ‘j’ to 

market ‘k’ 

Xyk Xjjk/XDk quantity transported as a fraction of total market demand 

yj 1 if warehouse is located at ‘j’,0 otherwise 

3.1.3 Mathematical Formulation 

Now cost minimization problem for the SSCWLP can be written as the mixed integer 
programming problem as follows 

FORMULATION SI VI. 1(A) 


MIN 


Cijk.Xijk ~b Z fi-yi 

i j k j 


( 0 ) 


Subject to — 


Xijk = 1 Vt, j, k 


( 1 ) 


i j k 


^ ' Xyk — ‘S’/V i 
j k 


( 2 ) 
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( 3 ) 


Xijk > dk\fk 

i j 

Xijk>(Ni,j\k 

zz xyk < capjyjij 

i k 

- y^j 

i k 

» = (0,l)Vy 

So the optimal integer solution to the problem is obtained by solving the above mixed 
integer formulation. 

MINIMIZE (0) subject to (1) to (7) 

The first part of the objective function denotes the cost of transportation of the 
commodity from the plant to the market via the warehouse; the second part is the fixed 
cost of locating the warehouse. 

Constraint 1 ensures that demand is met at each market point. 2 ensure that plants 
supply within their supplying capacities; 3 ensures that the demand at each point is met, 
whereas 4 is the non-negativity constraint. Constraint 5 ensures that the warehouses are 
routing the commodity to the markets from the plants within their capacity which is 
known; equation 6 is the linking constraint and ensures that warehouse is located only if 
the quantity shipped from the warehouse point is greater than zero. Equation 7 is the 
integer constraint on y\. 



(5) 

( 6 ) 

(7) 
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Now we relax the constraint 7 and so the problem now become 
FORMULATION SIVl.lflBl 


MINIMIZE (0) subject to (1) to (6) and (8) where 

yj > 0\/j ( 8 ) 

This can be referred as the weak relaxation as termed so in the literature ( Krarup and 
Purzan, (1983) [12]; Erlenkotter (1978) [14] as it is on the similar lines for weak 
relaxations of SPLP (Simple Plant Location Problem); and for CPLP (Capacitated Plant 
Location Problem) as given in Comuejols et.al. (1991) [8], 

In the next few variation of the model constraint (6) is relaxed and replaced by other 
linking constraints. 

FORMULATION S 1 VI 2(B) 

50 when we replace the constraint (6) and replace it by (9) as given below the problem 
becomes 

MINIMIZE (0) 

Subject to (1) to (5), (8) and (9) 

51 Xijk + M (1 — yi) > 0V j 

i k 

ss Xijk + Myj ^ 0 V/ (9) 

i k 

X£x, yi -My,<0V/ 

i k 

Equation (9) is another form of the linking constraint and ensures that the warehouse is 
located only when there is quantity shipped through the warehouse is greater than zero. 
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We call this constraint the Big M constraint where M is a constant having a very large 
value. 

FORMULATION SI VI .3(B) 


So when we replace the constraint (6) and replace it by (10) as given below the problem 

becomes 

MINIMIZE (0) 

Subject to (1) to (5), (8) and (10) 


y, xyk < yjdk\/j, k 


( 10 ) 


Equation (10) is another form of linking constraint and also ensures that the warehouse 
is located only if quantity shipped through the warehouse is greater than zero. In the 
literature the above relaxation is known as the strong relaxation. 


FORMULATION SI VI .4(B) 

So when we replace the constraint (6) and replace it by (1 1) as given below the problem 

becomes 

MINIMIZE (0) 

Subject to (1) to (5), (8) and (1 1) 

^ Xijk -Mil- yj) < dk\fj , k 

i 

y xijk — Myj < 0 V /, k ( H ) 

i 

y xijk + Myj > 0V j, k 
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Equation (1 1) is another form of linking constraint and also ensures that the warehouse 
is located only if quantity shipped through the warehouse is greater than zero. This 
constraint is another variation of the Big M constraint where M is a constant having a 
very large value. 


FORMULATION S 1 VI .5(B) 


So when we replace the constraint (6) and replace it by (10) and (12) as given below the 
problem becomes 
MINIMIZE (0) 

Subject to (1) to (5), (8) and (10) and (12) 

J]xijk< yjsNiJ (12) 

k 

Equation (10) and (12) are another form of linking constraints and also ensures that the 
warehouse is located only if quantity shipped through the warehouse is greater than 
zero. They relate the quantity shipped from a plant ‘i’ to the market ‘k’ through 
warehouse ‘j’ to the supply from the plant and the demand at the market. This is 
another form of a strong relaxation. 


FORMULATION SI VI. 6(B) 


So when we replace the constraint (6) and replace it by (1 1) and (13) as given below the 
problem becomes 
MINIMIZE (0) 

Subject to (1) to (5), (8) and (1 1) and (13) 
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(13) 


X Xi J' k -Mil- yj) < stf i, j 

k 

^ xuk - Myj < 0 V i, j 

k 

^ Xy* + Myj > 0 V i , j 

k 

Equation (11) and (13) are another form of linking constraint and also ensures that the 
warehouse is located only if quantity shipped through the warehouse is greater than 
zero. They relate the quantity shipped from a plant ‘i’ to the market ‘k’ through 
warehouse ‘j’ to the supply from the plant and the demand at the market. This is 
another form of a Big M relaxation. 

FORMULATION SlV1.7fBi 

So when we totally remove the linking constraint (6) as given below the problem 

becomes 

MINIMIZE (0) 

Subject to (1) to (5), (8) 

This is done to see the bounds given by the problem without the linking constraint and 
to see the strength of yj in the constraint (5). 

For all the above formulations SI VI. 2(B) to SI VI. 7(B) we have the problems giving 
the optimal integer solution just by replacing the constraint (8) with constraint (7) thus 
giving formulations SI VI. 2(A) to S1V1.7(A). At the stage of computational experience 
the integer formulation are used to verify that all the formulations give the same optimal 
results. 
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3.2 Problem formulation (STYLE 1 and Variation 2) 


In this section we will enlist all the formulation with Style 1 same as the style in the 
above section but only variation being in constraint (5). In this variation all the 
problems (from 1(A) to 6(B)) have the identical structures except that the constraint (5) 
is added in a different form as given below. 

^^Xijk < capjVj ► (14) 

i k 

So the integer formulation for the variation 2 is as follows: 

FORMULATION S Y V2. 1 (A) 

MINIMIZE (0) subject to (1) to (4) ; (6) to (7) and (14) 

Now we relax the constraint 7 and so the problem now becomes 

FORMULATION S1V2.UB) 

MINIMIZE (0) subject to (1) to (4), (6), (8) and (14) 

In the next few variation of the model constraint (6) is relaxed and replaced by other 
linking constraints. 

FORMULATION S1V2.2(B) 

So when we replace the constraint (6) and replace it by (9) as given below the problem 

becomes 

MINIMIZE (0) 

Subject to (1) to (4), (8), (9) and (14) 
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FORMULATION S1V2.3(B) 


So when we replace the constraint (6) and replace it by (10) as given below the problem 

becomes 

MINIMIZE (0) 

Subject to (1) to (4), (8), (10) and (14) 

In the literature the above relaxation is known as the strong relaxation. 

FORMULATION SlV2.4fBt 


So when we replace the constraint (6) and replace it by (1 1) as given below the problem 

becomes 

MINIMIZE (0) 

Subject to (1) to (4), (8), (1 1) and (14) 

FORMULATION S1V2.5(B) 

So when we replace the constraint (6) and replace it by (10) and (12) as given below the 
problem becomes 
MINIMIZE (0) 

Subject to (1) to (4), (8), (10), (12) and (14) 

FORMULATION S1V2.6IB) 

So when we replace the constraint (6) and replace it by (1 1) and (13) as given below the 
problem becomes 
MINIMIZE (0) 

Subject to (1) to (4), (8), (1 1), (13) and (14) 


For all the above formulations S1V2.2(B) to S1V2.6(B) we have the problems giving 
the optimal integer solution just by replacing the constraint (8) with constraint (7) thus 
giving formulations S1V2.2(A) to S1V2.6(A). At the stage of computational experience 
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the integer formulation are used to verify that all the formulations give the same optimal 
results. 


3.3Problem formulation (STYLE 2 and Variation 1) 

In this section we will enlist the various formulations of style 2 as adopted by Sharma 
(1991) [2]. Here the transportation of the commodity is done in two stages. In first stage 
the commodity is considered to be transferred form plant to the warehouse and then in 
next stage from the warehouse to the market. 

3.3.1 Constants Definition 

Dk Demand for the commodity at market ‘k’ 

dk Dk/EDk demand at market ‘k’ as a fraction of total market demand. 

Sj Supply available at plant ‘i’ of the commodity 

Sj Sj/EDk supply available at plant T as a fraction of the total market 

demand 

fj Fixed cost of locating a warehouse at ‘j’ 

Cpwij Unit cost of transporting goods from plant ‘i’ to warehouse ‘j 5 

Cwmjk Unit cost of transporting goods from warehouse : j’ to market ‘k’ 

CAPj Capacity of a warehouse ‘j ’ 

capj CAPj/EDk capacity of a warehouse at ‘j’ as a fraction of the total market 

demand 

3.3.2 Variable Definition 

Xpwij Quantity of commodity transported from plant ‘i’, to warehouse ‘j’. 

Xp W ij Xpwij/XDk quantity transported as a fraction of total market demand. 

Xwmjk Quantity of commodity transported from warehouse ‘j’ to market ‘k’. 

Xwmjk Xwmjk/XDk quantity transported as a fraction of total market demand, 

yj 1 if warehouse is located at ‘j’,0 otherwise 
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33.3 Mathematical Formulation 

Now cost minimization problem for the SSCWLP can be written as the mixed integer 


programming problem as follows 


FORMULATION S2V1 . 1 (A) 



Xpwij .Cpwij + / f Xwmjk .Cwmjk + 


MIN *' j j k 


SUBJECT TO 


^ i j Xpwij < Si\/ 1 

(15) 

j 

y ' Xwmjk ^ dkXf k 

(16) 

J 

Xpwij = 1 

i i 

(17) 

1 J 

j k 

(18) 

Xpwij 

(19) 

Xwmjk ^ 0V/ 5 k 

(20) 

'y ^ Xpwij ~ ^X\imjJc\fj 

i k 

(21) 

^ Xpwij < capjyj\fj 

i 

(22) 

P 

VI 

k 

• w- 

(23) 

»=(o>i)Vj 

(24) 


X fjyj 


( 00 ) 


So the optimal integer solution to the problem is obtained by solving the above mixed 
integer formulation. 



MINIMIZE (00) subject to (15) to (24) 


The first part of the objective function denotes the cost of transportation of the 
commodity from the plant to the warehouse the second part is the cost of transportation 
of the commodity from the warehouse to the market; and the third part is the fixed cost 
of locating the warehouse. 

Constraint 15 ensures that plants supply within their supplying capacities; 16 ensures 
that the demand at each point is met, whereas 17 and 1 8 ensure that the demand is met 
at each point. Constraints 19 and 20 are the non-negativity constraint. Constraint 21 is 
the flow balance constraint and ensures that the commodity coming into the warehouse 
is passed to the markets. Constraint 22 ensures that the warehouses are routing the 
commodity to the markets from the plants within their capacity which is known; 
equation 23 is the linking constraint and ensures that warehouse is located only if the 
quantity shipped from the warehouse point is greater than zero. Equation 24 is the 
integer constraint on yi. 

Now we relax the constraint 24 and so the problem now becomes 


FORMULATION S2V 1.1(B) 


MINIMIZE (00) subject to (15) to (23) and (25) where 

yj^ovj ( 25 ) 
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This can be referred as the weak relaxation as termed so in the literature ( Krarup and 
Purzan, (1983) [12]; Erlenkotter (1978) [14] as it is on the similar lines for weak 
relaxations of SPLP (Simple Plant Location Problem); and for CPLP (Capacitated Plant 
Location Problem) as given in Comuejols et.al. (1991) [8], 

In the next few formulations of the model constraint (23) is relaxed and replaced by 
other linking constraints. 

FORMULATION S2 VI .2(B) 

When we replace the constraint (23) and replace it by (26) as given below the problem 

becomes 

MINIMIZE (00) 

Subject to (15) to (22), (25) and (26) 

Xpwij + M(l- yi) > 0 V j 

i 

X Xpwij + Myj > 0 Vy (26) 

i 

yi Xpwij — Myj < 0 Vy 

i 

Equation (26) is another form of the linking constraint and ensures that the warehouse is 
located only when there is quantity shipped to the warehouse is greater than zero. We 
call this constraint the Big M constraint where M is a constant having a very large 
value. 

FORMULATION S2VL3(B) 

So when we replace the constraint (23) and replace it by (27) as given below the 
problem becomes 
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MINIMIZE (00) 

Subject to (15) to (22), (25) and (27) 

Xwmjk yjdkS/j, k ( 27 ) 

Equation (27) is another form of linking constraint and also ensures that the warehouse 
is located only if quantity shipped from the warehouse is greater than zero. It links the 
demand at market ‘k’ and the quantity to be transported. In the literature the above 
relaxation is known as the strong relaxation. 


FORMULATION S2VL4( 


So when we replace the constraint (23) and replace it by (28) as given below the 
problem becomes 
MINIMIZE (00) 

Subject to (15) to (22), (25) and (28) 


Xwmjk 1S/L (1 yk) ^ dkVj, k 

Xwmjk -Myj <0\/j,k 
Xwmjk + Myj > 0 V /, k 


Equation (28) is another form of linking constraint and also ensures that the warehouse 
is located only if quantity shipped from the warehouse is greater than zero. This 
formulation is another variation of the Big M constraint where M is a constant having a 
very large value. 
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FORMULATION S2VL5(B'i 


So when we replace the constraint (23) and replace it by (27) and (29) as given below 
the problem becomes 
MINIMIZE (00) 

Subject to (15) to (22), (25) and (27) and (29) 

Xpmj<yjSi\/i,j (29) 

Equation (27) and (29) are another form of linking constraint and also ensures that the 
warehouse is located only if quantity shipped through the warehouse is greater than 
zero. They relate the quantity shipped from a plant T to the market ‘k’ through 
warehouse ‘j’ to the supply from the plant and the demand at the market. This is 
another form of a strong relaxation. 

FORMULATION S2VL6(B) 

When we replace the constraint (23) and replace it by (28) and (30) as given below the 
problem becomes 
MINIMIZE (00) 

Subject to (15) to (22), (25) and (28) and (30) 

Xpwij Ad (1 yj) ^ si\f i, j 

Xpwij — Myj < 0 V z, j ( 30 ) 

Xpwij + Myj d. 0 

Equation (28) and (30) are another form of linking constraint and also ensures that the 
warehouse is located only if quantity shipped to and from the warehouse is greater than 
zero. They relate the quantity shipped from a plant T to the market ‘k’ through 
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warehouse £ j’ to the supply from the plant and the demand at the market. This is 
another form of a Big M relaxation. 


FORMULATION S2VL7£B} 

So when we totally remove the linking constraint (23) as given below the problem 

becomes 

MINIMIZE (00) 

Subject to (15) to (22), (25) 

This is done to see the bounds given by the problem without the linking constraint and 
to see the strength of yj in the constraint (8). 

For all the above formulations S2V 1.2(B) to S2V1.7(B) we have the problems giving 
the optimal integer solution just by replacing the constraint (25) with constraint (24) 
thus giving formulations S2V1.2(A) to S2V1.7(A). At the stage of computational 
experience the integer formulation are used to verify that all the formulations give the 
same optimal results. 

3.4 Problem formulation (STYLE 2 and Variation 2) 

In this section we will enlist all the formulation with Style 2 same as the style in the 
above section but only variation being in constraint (22). In this variation all the 
problems (from 1(A) to 6(B)) have the identical structures except that the constraint 
(22) is added in a different form as given below. 

Xpwij ^ CdJ?j\f J (31) 

i 
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So the integer formulation for the variation 2 is as follows: 

FORMULATION S2V2.UA) 

MINIMIZE (00) subject to (15) to (21), (23) to (24) and (31) 

Now we relax the constraint 24 and so the problem now becomes 

FORMULATION S2V2.1tB) 

MINIMIZE (00) subject to (15) to (21),(23), (25) and (31) where 

In the next few formulations of the model constraint (23) is relaxed and replaced by 
other linking constraints. 


FORMULATION S2V22IB) 

When we replace the constraint (23) and replace it by (26) as given below the problem 

becomes 

MINIMIZE (00) 

Subject to (15) to (21),(25), (26) and (31) 

FORMULATION S2V2.3fBt 

So when we replace the constraint (23) and replace it by (27) as given below the 
problem becomes 
MINIMIZE (00) 

Subject to (15) to (21), (25), (27) and (31) 

FORMULATION S2V2.4fB) 

So when we replace the constraint (23) and replace it by (28) as given below the 
problem becomes 
MINIMIZE (00) 

Subject to (15) to (21), (25), (28) and (31) 
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FORMULATION S2V2.S(B'i 

So when we replace the constraint (23) and replace it by (27) and (29) as given below 
the problem becomes 
MINIMIZE (00) 

Subject to (15) to (21), (25), (27), (29) and (31) 

FORMULATION S2V2.61B) 

When we replace the constraint (23) and replace it by (28) and (30) as given below the 
problem becomes 
MINIMIZE (00) 

Subject to (15) to (21), (25), (28), (30) and (31) 

For all the above formulations S2V2.2(B) to S2V2.6(B) we have the problems giving 
the optimal integer solution just by replacing the constraint (25) with constraint (24) 
thus giving formulations S2V2.2(A) to S2V2.6(A). At the stage of computational 
experience the integer formulation are used to verily that all the formulations give the 
same optimal results. 

Hence in this chapter we enlist two styles of formulations with their two variations and 
with all possible relaxations. 
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3.5 Summary 

In this section we give the summary of the formulations given in this chapter and 
give the nomenclature to be used in the subsequent sections: 


For both the Styles and Variations the formulations can be summarized as below 

*.1B: Weak relaxation of SSCWLP 
*.3B: SRS of one kind (yjdk>= £j x^k) 

*.5B: SRS of another kind (yjdk>= Si xp & yjSj>= I* Xjjk) 

*.2B: Big-M of first kind 

*.4B: Big-M of another kind 

*.6B: Two Big-M relaxations combined 

Variation 1: capjyj>=Ejk xp 

Variation 2: capj>=Ej k Xjjk 


42 



Chapter 4 

STRUCTURE OF EMPERICAL INVESTIGATION CARRIED 


4.1 Introduction 

We have tested the two styles of formulation first given by Geofffion and Graves (1974) 
[1] and the second one used by Sharma (1991) [2], For each of these styles we have two 
variations; the first one having a stronger relaxation for the capacity constraint and the 
other one having not so strong relaxation for the capacity constraint. For all these styles 
and variations there are seven different formulations in each category obtained by 
relaxing the linking constraint as listed down in chapter 3. 

We have tested the various relaxations for solving SSCWLP on a large variety of 
problems so as to see the effect of various parameters on the results. We have tried both 
small and large size problems for the warehouse location model and in this chapter we 
give the computational results as obtained by solving the problems. We have used the 
optimization package LINGO for solving the problems for various relaxations of 
SSCWLP. All the problems have been randomly generated and the results for the 
objective values and the number of iterations have been listed down. 

The problem testing phase has included the problems to be tested in the following 
manner. 
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4.2 Computational Details 


TABLE 4. 1 (PROBLEMS FOR PROBLEM SIZE 5X5X5) 



Style 

type 

Variation 

type 

Capacity 

Supply 

Relaxation 

IIPPI 

5X5X5 

Style 1 

Variation 1 

Ecapj= 2.0 

Espl,!. 5,2,5, 
10,20,30,40,50 

For each 
relaxation 

90 

5X5X5 

Style 1 

Variation2 

Eeap,-= 2.0 

Esj=T,1.5,2,5, 

10,20,30,40,50 

For each 
relaxation 


5X5X5 

Style 2 

Variation 1 

o 

<N 

II 

£ 

£ 

Es;=l,1.5,2,5, 

10,20,30,40,50 

For each 
relaxation 

90 

5X5X5 



Ecapj= 2.0 

Espl, 1.5 ,2,5, 
10,20,30,40,50 

For each 
relaxation 

90 

5X5X5 



Scapj= 6.0 

Espl, 1. 5,2,5, 
10,20,30,40,50 

For each 
relaxation 

90 

5X5X5 

B9I 


Ecapj= 6.0 

Espl, 1.5,2, 5, 
10,20,30,40,50 

For each 
relaxation 

90 

5X5X5 

Style 2 

Variationl 

o 

VO 

II 

£ 

Espl, 1.5 ,2,5, 
10,20,30,40,50 

For each 
relaxation 

90 

5X5X5 

Style 2 

Variation2 

Ecapj= 6.0 

EsrU.5,2,5, 

10,20,30,40,50 

For each 
relaxation 

90 

TOTAL NO. OF PROBLEMS 

720 


The Problem size 5X5X5 means that there are 5 plants, 5 warehouses to choose from 
and 5 demand sites. 

Similarly problems of size 10X10X10 and 20X20X20 are tried. 


TABLE 4.2(PROBLEMS FOR PROBLEM SIZE 10X10X10) 


Problem 

size 

Style 

type 

Variation 

type 

Capacity 

Supply 

Relaxation 

No of 

problems 

10X10X10 

Style 1 

Variationl 

Ecapj= 2.0 

Esi=l,1.5,2,5, 

10,20,30,40,50 

For each 
relaxation 

90 

10X10X10 

Style 1 

Variation2 

Ecapj= 2.0 

Esi=l,1.5,2,5, 

10,20,30,40,50 

For each 
relaxation 

90 

10X10X10 

Style 2 

Variationl 

Ecapj= 2.0 

Esj=l,1.5,2,5, 

10,20,30,40,50 

For each 
relaxation 

90 

10X10X10 

Style 2 

Variation2 

Ecapj= 2.0 

Esj=T,1.5,2,5, 

10,20,30,40,50 

For each 
relaxation 

90 

10X10X10 

Style 1 

Variationl 

Ecapj= 6.0 

Espl, 1.5, 2,5, 
10,20,30,40,50 

For each 
relaxation 

90 
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■ 

Style 1 

Variation2 

Ecapj= 6.0 

2sj=l,1.5,2,5, 

10,20,30,40,50 

For each 
relaxation 

90 

10X10X10 

Style 2 

Variationl 

o 

kd 

II 

f 

Esj=l,1.5,2,5,10 

,20,30,40,50 

For each 
relaxation 

90 





Sspl.l. 5,2,5, 
10,20,30,40,50 

For each 
relaxation 

90 

TOTAL N 

O. OF PI 

ROBLEMS 



TABLE 4.3(PROBLEMS FOR PROBLEM SIZE 20X20X20) 



Style 

type 

Variation 
Jype 

Capacity 

Supply 

Relaxation 

VSBS9 





£s;=l,l. 5,2,5, 
10,20,30,40,50 

For each 
relaxation 

90 

20X20X20 

Style 1 

Variation2 

O 

CN 

II 

§ 

8 

Xsj=l,1.5,2,5, 

10,20,30,40,50 

For each 
relaxation 

90 

20X20X20 

Style 2 

Variationl 

Ecapj= 2.0 


For each 
relaxation 

90 





Esj=l,l. 5,2,5, 
10,20,30,40,50 

For each 
relaxation 

90 

20X20X20 

Style 1 

Variationl 

Scapj= 6.0 

2sf=l,l. 5,2,5, 
10,20,30,40,50 

For each 
relaxation 

90 

20X20X20 

Style 1 

Variation2 

Xcapj= 6.0 

Ssj=l,l. 5,2,5, 
10,20,30,40,50 

For each 
relaxation 

90 

20X20X20 

Style 2 

Variationl 

Scapj= 6.0 

Esi=l,l. 5,2,5, 
10,20,30,40,50 

For each 
relaxation 

90 

20X20X20 

Style 2 

Variation2 

2capj= 6.0 

Ssi=l,l. 5,2,5, 
10,20,30,40,50 

For each 
relaxation 

90 

TOTAL N 

O. OF P] 

ROBLEMS 

720 


So in all 2160 different problems were tried to computationally establish the results and 
to have an effective analysis at later stages. 


Problems of bigger sizes could not be tried as the optimization package LINGO 
available in the labs is an educational package and has limitations on the number of 
variables and constraints in a problem that could be tried on the package. 
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4.3 Results 


In Appendix 1 the results obtained during the solving of various size problems ranging 
form 5X5X5 to 20X20X20 are given for various formulations. The objective value 
obtained by solving each of the problem and the number of iterations taken is shown. 


46 



Chapter 5 


ANALYSIS OF DATA AND RESULTS 


5.1 Introduction 

In this chapter statistical analysis of the data obtained in chapter 4 is done so as to bring 
forward some conclusive results from the work. The organization of the chapter is as 
follows. 

In section 5.2 we will list down the results for the main findings of the work; that the 
results obtained by Big M formulations (S1V1.2B, S1V1.4B, S1V1.6B) give 
statistically significant bounds obtained by other relaxations (S1V1.1B, S1V1.3B, 
S1V1.5B) for problems with Lcap,= 2.0 and for both the styles. Also in this section 
findings for case Ecapj= 6.0 is included. In section 5.3 we give the analysis results to 
show that the bounds obtained by S1V1.1B, S1V1.3B, S1V1.5B are same and have no 
statistically significant difference. In section 5.4 we will list down the results of analysis 
proving that the bounds obtained by Variation 1 for each problem category has better 
bounds than obtained by Variation 2. In section 5.5 we will give the results showing that 
in most circumstances Style 2 gives the same bounds as Style 1 with less number of 
iterations that is Style 2 is computationally less expensive than Style 1. Conclusions are 
listed down after each of the results establishing it. In the following sections all the T- 
tests and F-test have been done considering the null hypothesis that the difference 


between the means is zero. 



5.2.1 Statistical Significance of Big M Formulations (Case Ecapj= 2.0) 

This section deals with the most important finding of the work that when Ecapj = 2.0 

then for both the styles and for Variation 1 the Big M formulations (SIV1.2B, 
S1V1.4B,S1V1.6B,S2V1.2B, S2V1.4B, S2V1.6B) give the same bounds which are 
better than the bounds given by other relaxations (S1V1.1B, S1V1.3B, S1V1.5B, 
S2V1.1B, S2V1.3B, S2V1.5B). Since the bounds obtained by S1V1.1B, S1V1.3B, 
S1V1.5B, and S2V1.1B, S2V1.3B, S2V1.5B are same with no statistically significant 
difference we have carried out the T-test on S1V1.1B and S1V1.2B and also on 
S2V1.1B and S2V1.2B so as to establish that Bound (.2B) are significantly greater than 
Bound (.IB) 


Hence establishing that 


Bound (.2B, .4B, .6B) are significantly greater than Bound (.IB, ,3B, .5B) 


The results of the T-test is given in 
TABLE 5.1(T-test between .IB’s and .2B’s) 


Problem description 

T calculated 



5X5X5 SI VI .IB AND S1V1.2B 

-14.08 

5X5X5 S2V1 .1 B AND S2V1 .2B 

-14.10 

5X5X5 SI V2.1 B AND SI V2.2B 

77.67 

5X5X5 S2V2.1 B AND S2V2.2B 

. 77.66 





10X10X10 S1V1.1B AND S1V1.2B 

-16.59 

10X10X10 S2V1.1B AND S2V1.2B 

-16.59 

10X10X10 S1V2.1B AND SI V2.2B 

95.68 

10X10X10 S2V2.1B AND S2V2.2B 

95.74 





20X20X20 SI VI .IB AND SI VI .28 

-24.31 

20X20X20 S2V1 .1 B AND S2V1 ,2B 

-24.31 

20X20X20 S1V2.1B AND S1V2.2B 

125.63 

20X20X20 S2V2.1 B AND S2V2.2B 

125.64 
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The above values for T-calculated can be compared with the following values of T- 
critical at the different significance levels (for set of 90 data). 

Table 5.2 (T-critical values data size 90) 


Significance 

level 

Alpha= 0.05 

Alpha= 0.01 

Alpha= 0.005 

T-critical(one 

tail) 

1.662156 

2.368979 

2.632205 

T-criticaI(two 

tailed) 

1.986978 

2.632205 



To support the fact that the Big M formulations give significantly better bounds is 
strengthened by the table 5.3 given below indicating the % improvement in the bounds 
given by Big M formulations over other relaxations. 

Table 5.3 (% Improvement of Big-M formulation over SRS formulation of 
SSCWLP in Variation 1) 


Problem Description 

% Improvement 

5X5X5 Ecapj= 2.0 SI VI. IB AND S1V1.2B 

10.21% 

5X5X5 Ecapj= 2.0 S2V1.1B AND S2V1.2B 

10.24% 

10X10X10 Ecap,= 2.0 SI VI. IB AND S1V1.2B 

24.12% 

10X10X10 Ecap,-= 2.0 S2V1.1B AND S2V1.2B | 

24.18%. 

20X20X20 Ecap,= 2.0 SI VI. IB AND S1V1.2B 


20X20X20 Ecapj= 2.0 S2V1.1B AND S2V1.2B 

38.84% 


From the data obtained above in the analyses phase we can infer the following from the 
tables 5.1, 5.2, 5.3: 

1 . for Ecapj= 2.0 the bounds obtained in Variation 1 and both the Styles for 
formulations .2B, .4B, .6B are significantly better than those obtained by other 
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relaxations. Big-M constraint really boosts the bounds in Variation 1 
(capjyj>=SikXiji) for case leap j- 2.0. 

2. The improvement increases with the problem size so we can hope for more 
significant improvements for bigger problem sizes which tend to the real world 
problems. 

3. The bounds obtained by in Variation 2 and both the Styles for formulations ,2B, 
.4B, and .6B are worse than bounds obtained by formulations .IB, .3B, and .5B. 
Hence we see that the strength of yj in the capacity constraint for both the Styles 
and all the formulations is significant and so in Variation 2 of all the 
formulations we get very poor results for formulations .2B, .4B, and .6B. Big-M 
is worse than SRS & WRS in Variation 2 (capj>= HikXijk) for £capj= 2.0. 

Conclusion 1: 

For both the Styles and both Variations we get same bounds for .2B, .4B and .6B. 

Conclusion 2: 

For both the styles and Variation 1 with low capacities (Zcapj= 2.0) and all problem 

sizes 

Bound (.2B, ,4B, .6B) are significantly greater than Bound (.IB, .3B, .5B) 

5.2.2 Statistical Significance of Big M Formulations (Case Scapj= 6.0) 

In this section the results for the T-test carried for case Ecapj= 6.0 between the .IB 

formulations and ,2B for both the Styles and Variations yield the results is shown in 

Table 5.4 
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TABLE 5.4(T-test between .IB’s and .2B’s) 


Problem description 

T calculated 



5X5X5 S1V1.1B AND SI VI .28 

24.23 

5X5X5 S2V1.1B AND S2V1.2B 

24.22 

5X5X5 S1V2.1B AND S1V2.2B 

64.48 

5X5X5 S2V2.1B AND S2V2.2B 

64.90 





10X10X10 S1V1.1B AND S1V1.2B 

17.70 

1 0X10X1 0 S2V1 . 1 B AND S2V1 .2B 

17.76 

10X10X10 S1V2.1B AND S1V2.2B 

89.81 

10X10X10 S2V2.1B AND S2V2.2B 

89.78 





20X20X20 SI VI. IB AND S1V1.2B 

0.47 

20X20X20 S2V1 .1 B AND S2V1 .2B 

0.26 

20X20X20 S1V2.1B AND S1V2.2B 

141.79 

20X20X20 S2V2.1B AND S2V2.2B 

141.56 


From the above table we can enlist the above inferences: 


1. For Ecapj= 6.0 the bounds obtained by in Variation 2 and both the Styles for 
formulations ,2B, ,4B, and .6B are worse than bounds obtained by 
formulations .IB, ,3B, and ,5B. This implies that Big-M is ineffective in 
boosting bounds in Ecapj= 6. 0. 

2. For Ecapj= 6.0 the bounds obtained in Variation 1 and both the styles for 
formulations .2B, .4B, .6B are worse than bounds obtained by formulations 
.IB, .3B, .5B in case when the problem size is 5X5X5 and 10X10X10. But 
in problem size 20X20X20 there is no significant difference between bounds 
obtained by the two formulations. We can say that with Ecapj= 6.0 the 
problems for small problem sizes can be treated as Uncapacitated problem 
where any one warehouse location is sufficient to fulfill all the demands. 
Hence the bounds obtained by the formulations .IB, .3B, .5B perform better 


** O S i ^ 

& ***** 
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as yj takes value 1 for the LP relaxation. Big-M fails to boost the bounds that 
is contradictory to the caseof what happened in Ecap?= 2.0. 

Hence with the help of the data in chapter 4 and above tables we can give following 
conclusions 


Conclusion 3: 

For both the styles and both variations we get same bounds for .26, .46 and .66. 
Conclusion 4: 

For both the styles and variation 1 with Zcapj= 6.0 (high capacities) 

Hounds (.16, .36, .56) are significantly greater than 6ounds (.26, .46, .66) for small 
problem sizes (5X5X5 and 10X10X10) 

And 

Hounds (.16, .36, .56) no significant difference Hounds (.26, .46, .66) for larger 
problem sizes (20X20X20) 


5.3 Statistically insignificant difference of Bounds for .IB, .3B and .5B 

This section shows the analysis of the bounds obtained by the weak and the two strong 
relaxations as they are referred in chapter three. Though there is some improvement in 
the bounds given by .3B and .56 over the bound given by .16 for both the Styles and 
Variations and both capacity levels of Ecapj= 2.0 and Ecapj= 6.0. But when an F- test is 
done on the data it is seen that the difference is statistically insignificant. The results of 


the F-test done is shown as below in Table 5.5 



Table 5.5(F-Test between .IB, 3B and .5B) 


Problem Description I F -calculated 


5X5X5 Scapr 2.0 S1V1.1B,.3B,.5B 

8.63E-11 

5X5X5 £capj= 2.0 S2V1.1B..3B..5B 

0.000113 

5X5X5 Scapr 2.0 S1V2.1B,.3B,.5B 

0.006605 

5X5X5 Scapr 2.0 S2V2.1B,.3B,.5B 

0.00855 




0.000532 

0.000448 

9.24E-05 

0.170664 



10X10X10 Ecapj= 2.0 S1V1.1B..3B..5B 

-4.70E-13 

10X10X10 Scapj= 2.0 S2V1 .1 B..3B..5B 

3.59E-05 

10X10X10 2capj= 2.0 S1V2.1B..3B..5B 

0.012838 

10X10X10 Scapj= 2.0 S2V2.1B..3B..5B 

0.016373 



10X10X10 £capj= 6.0 S1V1.1B..3B..5B 8.20E-08 

10X10X10 Scapj= 6.0 S2V1.1B,.3B,.5B 0.000782 

10X10X10 Scapr 6.0 S1V2.1B,.3B,.5B 0.00335 

10X10X10 Zcapj= 6.0 S2V2.1B,.3B,-5B | 0.003338 



20X20X20 Ecapj= 2.0 S1V1.1B,.3B,.5B 

2.80E-05 

20X20X20 Scapj= 2.0 S2V1.1B..3B..5B 

7.40E-06 

20X20X20 Ecapj= 2.0 S1V2.1B,'.3B,.5B 

0.04989 

20X20X20 Zcapp 2.0 S2V2.1B,.3B,.5B j 

0.052704 


20X20X20 Scapr 6.0 S1V1.1B..3B..5B 2.56E-07 


20X20X20 Scapr 6.0 S2V1.1B,.3B,.5B 

0.000484 

20X20X20 Zcapj= 6.0 S1V2.1B..3B..5B 

0.017559 

20X20X20 Scapr 6.0 S2V2.1B,.3B,.5B 

0.046795 









































The above values for F-calculated can be compared with the following values of F- 
critical at the different significance levels (for set of 90 data). 

Table 5.6 (F-critical values) 


Significance 

level 

Alpha= 0.05 

Alpha= 0.01 

Alpha= 0.005 

F-critical 

3.029598 

4.685546 

5.4049 


For the data given in tables 5.5 and 5.6 we can infer that there is no significant 
difference between the bounds obtained by the formulations .IB, .3B and .5B for both 
the styles and both the variations. 

Hence with the help of the data in chapter 4 and above tables we can give following 

conclusions 

Conclusion 5: 

There is no significant difference between Bound (.IB), Bound (-3B) and Bound (.5B) 
for both Styles and both the Variations. 

5.4 Better Bounds by Variation 1 than Variation 2 

This section shows the analysis of the bounds obtained by the all the relaxations of 
Variation 1 and Variation 2 and both capacity levels of Ecapj= 2.0 and Scapj= 6.0. It is 
seen from the analysis that bounds obtained by Variation 1 is significantly better that 
bounds obtained by relaxations in Variation 2 for both the Styles. This can be attributed 
to the fact that the addition of yj in capacity constraint gives a stronger relaxation for all 
formulations than the formulations without y, in capacity constraint. A T-test was done 
on data set one by one for the corresponding formulation in Variation 1 and Variation 2 
for both the capacity levels. Tables 5.7, 5.8, 5.9 given below establish the above fact. 
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Table 5.7(T-Test between Bounds of Variation 1 and Variation 2 for All 
Relaxations, Both Styles and both Capacity Levels; Problem Size 5X5X5) 


Problem Description 

T-Calcu!ated 



5X5X5 2capi= 2.0 S1V1.1B AND S1V2.1B 

18.31 

5X5X5 Scapf 2.0 S1V1.2B AND S1V2.2B 

50.83 

5X5X5 Sca Pi = 2.0 S1V1.3B AND S1V2.33 

18.27 

5X5X5 Scapr 2.0 S1V1.4B AND S1V2.43 

50.83 

5X5X5 £ca Pi = 2.0 S1V1.5B AND S1V2.5B 

18.27 

5X5X5 £capf= 2.0 SI VI .6B AND SI V2.6B 




5X5X5 £cap,= 2.0 S2V1 .1 B AND S2V2.1 B 

18.25 

5X5X5 Ecapr 2.0 S2V1 ,2B AND S2V2.2B 

50.83 

5X5X5 £cap,= 2.0 S2V1 ,3B AND S2V2.3B 

18.21 

5X5X5 Ecapr 2.0 S2V1 ,4B AND S2V2.4B 

50.83 

5X5X5 £ca Pi = 2.0 S2V1 ,5B AND S2V2.5B 

18.21 

5X5X5 £cap,= 2.0 S2V1 .6B AND S2V2.6B 

50.83 



5X5X5 £ca Pi = 6.0 S1V1.1B AND S1V2.1B 

5.50 

5X5X5 Eca Pf = 6.0 SI VI .2B AND SI V2.2B 

41.82 

5X5X5 Scapf* 6.0 S1V1.3B AND S1V2.3B 

5.53 

5X5X5 Ecapf* 6.0 SI VI .4B AND SI V2.4B 

47.90 

5X5X5 Ecapr 6.0 S1V1.5B AND S1V2.5B 

5.40 

5X5X5 2capi= 6.0 S1V1.6B AND S1V2.6B 



■■■■■■■ 

5X5X5 Ecap,= 6.0 S2V1 .1 B AND S2V2.1B 

5.40 

5X5X5 Ecapr 6.0 S2V1 .2B AND S2V2.2B 

41.89 

5X5X5 £cap,= 6.0 S2V1 ,3B AND S2V2.3B 


5X5X5 Xcapr 6.0 S2V1 .4B AND S2V2.4B 

■m 

5X5X5 Ecapr 6.0 S2V1 .5B AND S2V2.5B 


5X5X5 Ecapr 6.0 S2V1.6B AND S2V2.6B 



Table 5.8(T-Test between Bounds of Variation 1 and Variation 2 for AH 
Relaxations, Both Styles and both Capacity Levels; Problem Size 10X10X10) 


Problem Description 

T-Calculated 



10X10X10 Xcapi= 2.0 S1V1.1B AND S1V2.1B 

26.38 

10X10X10 Zcap,-= 2.0 S1V1.2B AND S1V2.2B 

52.88 

10X10X10 £ca P i= 2.0 S1V1.3B AND S1V2.3B 

26.34 

10X10X10 £ca P i= 2.0 S1V1.4B AND S1V2.4B 

52.88 

10X10X10 Ecap,— 2.0 SI VI .5B AND S1V2.5B 

26.34 

10X10X10 Ecapr 2.0 S1V1.6B AND S1V2.6B 

52.88 
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Problem Description 
10X10X10 Scapp 2.0 < 


10X10X10 Scapp 2.0 < 


S2V1.1B AND S2V2.1B 


S2V1.2B AND S2V2.2B 


10X10X10 Scapp 6.0 


10X10X10 Scapp 6.0 


10X10X10 Scapp 6.0 


10X10X10 Scapp 6.0 


10X10X10 Scapp 6.0 


10X10X10 Sea 


S2V1.1B AND S2V2.1B 


S2V1 .2B AND S2V2.2B 


S2V1 .3B AND S2V2.3B 


S2V1.4B AND S2V2.4B 


S2V1 .5B AND S2V2.5B 


S2V1 .6B AND S2V2.6B 


T-Ca leu fated 


10X10X10 Eca Pi = 2.0 S2V1 ,3B AND S2V2.3B 

26.34 

1 0X1 0X1 0 2cap,= 2.0 S2V1 .4B AND S2V2.4B 

52.91 

10X10X10 Scapp 2.0 S2V1 .5B AND S2V2.5B 

26.33 

10X10X10 Scapp 2.0 S2V1.6B AND S2V2.6B 

52.88 


10X10X10 Scapp 6.0 S1V1.1B AND S1V2.1B 

9.14 

10X10X10 Scapp 6.0 S1V1.2B AND S1V2.2B 

33.08 

10X10X10 Scapp 6.0 S1V1 ,3B AND S1V2.3B 

9.10 

10X10X10 Scapp 6.0 S1V1.4B AND S1V2.4B 

33.08 

10X10X10 Scapp 6.0 S1V1.5B AND S1V2.5B 

9.09 

10X10X10 Scapp 6.0 S1V1.6B AND S1V2.6B 

33.08 


Table 5.9(T-Test between Bounds of Variation 1 and Variation 2 for All 
Relaxations, Both Styles and both Capacity Levels; Problem Size 20X20X20) 



Problem Description 


T-Cafculated 


20X20X20 Sca Pi 
20X20X20 Sca Pi 
20X20X20 Xcapj 
20X20X20 Scap, 
20X20X20 Scap, 
20X20X20 Scan 


20X20X20 Scapp 2.0 
20X20X20 Scapp 2.0 
20X20X20 Scapp 2.0 
20X20X20 Scapp 2.0 
20X20X20 Sca Pi p 2.0 
20X20X20 Scapp 2.0 


SI VI. IB AND S1V2.1B 
S1V1.2B AND S1V2.2B 
S1V1.3B AND S1V2.3B 
S1V1.4B AND S1V2.4B 
S1V1.5B AND S1V2.5B 
S1V1.6B AND S1V2.6B 

S2V1.1B AND S2V2.1B 
S2V1.2B AND S2V2.2B 
S2V1-3B AND S2V2.3B 
S2V1.4B AND S2V2.4B 
S2V1.5B AND S2V2.5B 
S2V1.6B AND S2V2.6B 


20X20X20 Scapp 6.0 SI VI .1 B AND SI V2.1 B , 


20X20X20 Scapp 6.0 SI VI .2B AND SI V2.2B 
20X20X20 Scapp 6.0 S1V1.3B AND S1V2.3B 
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Problem Description 

1— Id" 1 ! H"?''" r| ""n 



20X20X20 Eca P j= 6.0 S 1 VI .4B AND S 1 V2.4B 

45.78 

20X20X20 2ca Pi = 6.0 S1V1.5B AND S1V2.5B 

12.28 

20X20X20 Xca Pi = 6.0 SI VI ,6B AND SI V2.6B 

45.78 



20X20X20 Scapp 6.0 S2V1.1B AND S2V2.1 B 

12.35 

20X20X20 Sca Pi = 6.0 S2V1 .2B AND S2V2.2B 

46.29 

20X20X20 2capf= 6.0 S2V1 .3B AND S2V2.3B 

12.27 

20X20X20 Ecapp 6.0 S2V1 .4B AND S2V2.4B 

46.29 

20X20X20 Ecapp 6.0 S2V1 .5B AND S2V2.5B 

12.28 

20X20X20 Ecapp 6.0 S2V1 .6B AND S2V2.6B 

46.29 


When we compare the calculated values for the T-statistic and the T-critical given in 
table 5.2 at the different significance levels (for set of 90 data) we can very well say that 
the bounds obtained by different formulations for Variation 1 is statistically significant 
than the bounds obtained by different formulation for Variation 2 for both the Styles and 
both capacity levels. 

Hence with the help of the data in chapter 4 and above tables we can give following 
conclusion. 


Conclusion 6: 

Bounds (S1V1.1B, ,2B, .3B, ,4B, .5B, .6B) 
Are significantly greater than 
Bounds (S1V2.1B, .2B, .3B, .4B, .5B, ,6B) 
And 

Bounds (S2V1.1B, .2B, .3B, .4B, .5B, ,6B) 
Are significantly greater than 
Bounds (S2V2.1B, .2B, .3B, .4B, .5B, ,6B) 


With corresponding comparison of different formulations. 
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5.5 Comparison of Number of Iterations for Both Styles 


This section shows the analysis of the number of iterations needed by the two Styles for 
all the relaxations of Variation 1 and Variation 2 and both capacity levels. It is seen 
from the analysis that the number of iterations required for a formulation in Style 2 is 
often less than number of iterations required for the corresponding formulation in Style 
1 . To establish this fact a T-test was done on the corresponding formulations for both 
the Variations and capacity levels. This can be attributed to the more complex nature of 
Style 1 and also to more no of variables and constraints compared to Style 2. The Tables 
5.10, 5. 11, 5. 12 given below illustrate the above fact. 

Table 5.10(T-Test on number of iterations for All corresponding Relaxations of 
Stylel and Style 2 and both Capacity Levels; Problem Size 5X5X5) 


Problem Description 

T-CaicuSated 



5X5X5 Xcap,= 2.0 S1V1.1A AND S2V1.1A 

-3.04 

5X5X5 2ca Pi = 2.0 SI VI .1 B AND S2V1 .1 B 

4.01 

5X5X5 Xca Pi = 2.0 SI VI .2B AND S2V1 .2B 

-0.85 

5X5X5 £cap,= 2.0 SI VI ,3B AND S2V1 .3B 

6.75 

5X5X5 Eca Pj = 2.0 SI VI .4B AND S2V1 .4B 

-5.03 

5X5X5 Eca Pi = 2.0 SI VI .5B AND S2V1 .5B 

5.94 

5X5X5 Ecap,= 2.0 SI VI .6B AND S2V1 .6B 

-4.91 

5X5X5 Ecap,= 2.0 S1V1.7B AND S2V1.7B 

7.65 



5X5X5 Eca P{ = 2.0 SI V2.1 A AND S2V2.1A 

-0.74 

5X5X5 Ecap;= 2.0 SI V2.1 B AND S2V2.1 B 

6.27 

5X5X5 Soap,— 2.0 SI V2.2B AND S2V2.2B 

4.67 

5X5X5 Sca Pi = 2.0 SI V2.3B AND S2V2.3B 

8.41 

5X5X5 2ca P i= 2.0 SI V2.4B AND S2V2.4B 

5.89 

5X5X5 Leap,- 2.0 SI V2.5B AND S2V2.5B 

7.62 

5X5X5 Scapr 2.0 SI V2.6B AND S2V2.6B 

6.91 



5X5X5 Ecap 

= 6.0 SI VI. 1 A AND S2V1.1A 

-5.15 

5X5X5 Ecap 

= 6.0 S1V1.tBANDS2V1.1B 

3.68 

5X5X5 Ecap 

= 6.0 S1V1.2B AND S2V1.2B 

2.51 

5X5X5 Ecap 

= 6.0 SI VI .3B AND S2V1 .38 

9.07 

5X5X5 Ecap 

= 6.0 S1V1.4B AND S2V1.4B 

-5.48 

5X5X5 Ecap 

= 6.0 SI VI .5B AND S2V1 .5B 

5 97 

5X5X5 Ecap 

= 6.0 S1V1.6BANDS2V1.6B 

-8.56 

5X5X5 Ecap 

= 6.0 S1V1.7B AND S2V1.7B 

3.89 
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Table 5.11(T-Test on number of Iterations for All corresponding Relaxations of 
Stylel and Style 2 and both Capacity Levels; Problem Size 10X10X10) 



Problem Description 


T-Calculated 


10X10X10 Scap,- 2.0 

SI VI .1 A AND S2V1 .1 A 

10X10X10 2cap;= 2.0 

S1V1.1B AND S2V1.1B 

10X10X10 £ca Pi = 2.0 

S1V1.2B AND S2V1.2B 

10X10X10 £cap,= 2.0 

S1V1.3BANDS2V1.3B 

10X10X10 Scapj=2.0 

S1V1.4B AND S2V1.4B 

10X10X10 Zcapj= 2.0 

S1V1.5B AND S2V1.5B 

10X10X10 2capj=2.0 

S1V1.6B AND S2V1.6B 

10X10X10 Sca Pi = 2.0 

S1V1.7B AND S2V1.7B 


1 0X1 0X1 0 Ecap,-= 2.0 

SI V2.1 A AND S2V2.1A 

10X10X10 Xca Pi = 2.0 

S1V2.1B AND S2V2.1B 

10X10X10 2capi= 2.0 

S1V2.2B AND S2V2.2B 

10X10X10 £capi=2.0 

SI V2.3B AND S2V2.3B 

10X10X10 Scapp 2.0 

S1V2.4B AND S2V2.4B 

10X10X10 2cap,-= 2.0 

S1V2.5B AND S2V2.5B 

10X10X10 2capj= 2.0 

S1V2.6B AND S2V2.6B 


10X10X10 Xca Pi = 6.0 


10X10X10 Scapp 6.0 


10X10X10 Scap;= 6.0 


10X10X10 Xcapr 6.0 


10X10X10 Xcapr 6.0 


10X10X10 Scapr 6.0 
10X10X10 Ecap,= 6.0 
10X10X10 Soap,- 6.0 


10X10X10 Xcapp 6.0 


10X10X10 Xcap,= 6.0 


10X10X10 Xcapr 6.0 


10X10X10 Zcapp 6.0 


1 0X1 0X1 0 Zcap,~ 6.0 
10X10X10 Eca P ;= 6.0 


10X10X10 Xcap,-= 6.0 


.1A AND 


.IB AND 


.2B AND 


,3B AND 


-4B AND 


,5B AND 
,6B AND 
.7B AND 


S2V1.1A 


S2V1.1B 


S2V1.2B 


S2V1.3B 


S2V1.4B 


S2V1.5B 

S2V1.6B 

S2V1.7B 


SI V2.1 A AND S2V2.1A 


S1V2.1B AND S2V2.1B 


S1V2.2B AND S2V2.2B 


S1V2.3B AND S2V2.3B 


S1V2.4B AND S2V2.4B 
S1V2.5B AND S2V2.5B 


S1V2.6B AND S2V2.6B 
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Table 5.12(T-Test on number of iterations for Ali corresponding Relaxations of 
Stylel and Style 2 and both Capacity Levels; Problem Size 20X20X20) 


Problem Description 

T-Calculated 



20X20X20 2capj= 2.0 SI VI. 1 A AND S2V1.1 A 

13.52 

20X20X20 Zcapr 2.0 SI VI .1 B AND S2V1 .1 B 

11.91 

20X20X20 Zcapi= 2.0 S1V1.2B AND S2V1.2B 

18.79 

20X20X20 £ca Pi = 2.0 SI VI ,3B AND S2V1 ,3B 

14.36 

20X20X20 Eca P j= 2.0 S1V1.4B AND S2V1.4B 

3.59 

20X20X20 Eca Pi = 2.0 SI VI ,5B AND S2V1 .58 

14.59 

20X20X20 2ca Pi = 2.0 S1V1.6B AND S2V1.6B 

4.41 

20X20X20 2ca Pj = 2.0 SI VI ,7B AND S2V1 ,7B 

8.80 



20X20X20 Scapr 2.0 SI V2.1 A AND S2V2.1A 

2.71 

20X20X20 Ecapp 2.0 SI V2.1 B AND S2V2.1 B 

11.75 

20X20X20 Zca P j= 2.0 SI V2.2B AND S2V2.2B 

20.16 

20X20X20 Ecapp 2.0 SI V2.3B AND S2V2.3B 

13.21 

20X20X20 Ecapp 2.0 SI V2.4B AND S2V2.4B 

8.71 

20X20X20 Ecapp 2.0 SI V2.5B AND S2V2.5B 

12.44 

20X20X20 Ecapp 2.0 SI V2.6B AND S2V2.6B 

13.35 



20X20X20 2cap,-= 6.0 S1V1.1A AND S2V1.1A 

12.23 

20X20X20 Ecapp 6.0 SI VI. IB AND S2V1.1B 

10.19 

20X20X20 2ca Pi = 6.0 SI VI ,2B AND S2V1 .2B 

18.77 

20X20X20 2cap,-= 6.0 SI VI .3B AND S2V1 .3B 

16.65 

20X20X20 Sca Pi = 6.0 SI VI .4B AND S2V1 ,4B 

2.94 

20X20X20 2ca Pi = 6.0 SI VI ,5B AND S2V1 ,5B 

7.49 

20X20X20 2cap,= 6.0 SI VI ,6B AND S2V1 .6B 

2.39 

20X20X20 2ca Pi = 6.0 SI VI ,7B AND S2V1 .7B 

7.16 



20X20X20 Zcap,-= 6.0 SI V2.1 A AND S2V2.1A 

7.17 

20X20X20 2cap,-= 6.0 S1V2.1B AND S2V2.1B 

6.99 

20X20X20 2cap ; = 6.0 SI V2.2B AND S2V2.2B 

16.95 

20X20X20 2cap,= 6.0 SI V2.3B AND S2V2.3B 

13.54 

20X20X20 2cap;= 6.0 SI V2.4B AND S2V2.4B 

10.00 

20X20X20 Scap,= 6.0 SI V2.5B AND S2V2.5B 

11.00 

20X20X20 2capj= 6.0 SI V2.6B AND S2V2.6B 

14 8" 


The above values for T-calculated can be compared with the following values of T- 
critical at the different significance levels (for set of 90 data) given in Table 5.2. 
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From the data obtained above in the analyses phase we can infer the following from the 
tables 5.10, 5.11, 5.12: 

1. For problem Size 20X20X20 for all formulations and both the variations the 
number of iterations for Style 1 is significantly higher than for Style 2. 

2. For problem size 5X5X5 and 10X10X10 there are some instances when the 
number of iterations for Style 2 is greater than Style 1, but in most formulations 
number of iterations for Style 1 is greater than Style2. The first case is 
encountered in formulations .4B and ,6B on some instances. 

3. For higher problem size the number of instances where number of iterations for 
Style 2 is greater than Style 1 is very low or almost absent and so it can be said 
that for higher problem sizes number of iterations for Style 1 is greater than 
Style2 for all the formulations. 

Hence with the help of the data in chapter 4 and above tables we can give following 
conclusion. 

Conclusion 7: 

Numbers of Iterations (Style 1) 

Are significantly greater than 
Number of Iterations (Style 2) 

We found that all seven conclusions are fully supported by the statistical tools. 



Chapter 6 

NEW FORMULATION 


Encouraged by the results obtained in chapter four for Big M formulations in this 
chapter we have tried a new formulation only for Style 1 and Variation 2 so as to see the 
effect of 

the bounds given when the Big M formulations are combined with the strong relaxation. 
So the new formulation can be listed as below: 

FORMULATION NEW 


So when we replace the constraint (6) and replace it by (10) and (1 1) as given below the 
problem becomes 

MINIMIZE (0) 

Subject to (1) to (4), (8), (10), (11) and (14) 

This was done to see if there is some further improvement in the bounds given by 
Formulation S1V2.3B after adding the Big M constraint with the constraint 10. 

For the analyses fifty problems for problem size 20X20X20 that are randomly generated 
are tried and results are given in Appendix 2. 

On the analyses of the result we see that although there is no improvement in the bounds 
given by the new formulation over the bounds given by S1V2.1B, S1V2.3B and 
S1V2.5B, but the number of iterations needed for the new formulation is significantly 
less than the number of iterations needed for the formulation S1V2.3B. Since in chapter 
5 we have already concluded that there is no significant difference between the bounds 
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given by .IB, .3B, .5B; we have done a T-test of the bounds given by 3B with bounds 
given by new formulation only and assume it to be applicable to .IB and .5B also. 

Thus we establish that there is no significant difference in the bounds given by the two 
formulations. 

The following table shows the results of the T-test done on the new formulation and 
SI V2.3B for both the objective value and the number of iterations. 

Table 6.1 (T-test for the New Formulation) 


OBJ VALUE 

T- Calculated 



20X20X20 CAP2 S1V2.3B AND NEW 

-0.16685 





NO OF ITERATIONS 




20X20X20 CAP2 S1V2.3B AND NEW 

3.227632 




If we compare the above T- calculated with the T-critical values at different significance 
levels we see that the Bounds given by the two formulations are not statistically 
different but the number of iterations needed by the new formulation Is significantly less 
than the number of iterations needed by SI V2.3B. 

The T- critical values at different significance levels for a data size of fifty is given as 
below: 

Table 6.2(T-criticaI data size 50) 


Significance 

level 

Alpha= 0.05 

Alpha= 0.01 

Alpha= 0.005 

T-critical(one 

tail) 

1.676551 

2.404886 

2.679953 

T-criticaI(two 

tailed) 

2.009574 

2.679953 

2.939742 
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Chapter 7 

DIRECTIONS FOR FUTURE RESEARCH 

Following future research topics are suggested: 

1. To develop a heuristic to establish bounds given by formulations S1V1.2B, .4B, 
.6B and S2V1.2B, .4B, .6B (Big M formulations) for low capacities as they give 
stronger relaxation bounds than the SRS. 

2. To give the theoretical proofs to establish the empirical finding in the present 
work. 
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STYLE 2, VARIATION 1(10X10X10. Zcap=6 n 


< 

<r* 

CD 

CD 

04 

to 

m 

m 

m 

z 

o 

*5 

3 

Z 

o 

H 

3 

3 

z 

o 

h- 

3 

3 

z 

o 

p 

3 

z 

o 

p 

3 

u» 

Z 

o 

J— 

3 

to 

Z 

o 

b- 

3 

2 

DC 

s 

cc 

s 

C£ 

w 

s 

3 

s 

I 

3 

s 

O 

UL 

o 

LL 

O 

LL. 

o 

LL 

DC 

o 

LL 

sr 

o 

LL 

DC 

o 

u_ 





nraSHUgai 


Esi=1.0 



FORMULATION7B 




















ON/S 
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STYLE 2, VARIATION 1(20X 20X20. 2cap=2.0 
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STYLE 1. VARIATION 2(20X20X20, £cap=6.0) 


S/NO 

FORMULATION 1A 

FORMULATION IB 

FORMULATION 2B 

FORMULATION 3B 

FORMULATION 4B 

FORMULATION 5B 

FORMULATION 6B 

I i 

iaaiiii3a 


M&Z MHga 




mssmuam 


Zsi=1.0 




8 
































































































^MMgawaia 


Esi=1.0 


KEM«H:1 


QKnit^jii 









































APPENDIX 3 



STYLE 1 VARIATION 1 


FORMULATION S1V1.1(A) 


SETS: 

SUPPLY/1. .20/ : S ; 

WAREHOUSE/1. .20/: CAP, F, Y; 

DEMAND/1. .20/ : D ; 

TRANSPORT (SUPPLY , DEMAND , WAREHOUSE ) : C , X ; 
DUMMY 1 (SUPPLY, DEMAND) ; 

DUMMY2 (WAREHOUSE, DEMAND) ; 

DUMMY3 (SUPPLY, WAREHOUSE) ; 

ENDSETS 


DATA: 


ENDDATA 

MIN=@SUM (TRANSPORT ( I, J,K) :C(I, J,K) *X(I, J,K) ) + 
@SUM (WAREHOUSE ( J) : F ( J) *Y { J) ) ; 


! SUBJECT TO 

@SUM (TRANSPORT ( I, J,K) :X(I, J, K) )=1; 

@FOR (SUPPLY (I) : @SUM (DUMMY2 (J,K) :X (I, J, K) ) <=S (I)); 

@FOR (DEMAND (K) :@SUM(DUMMY3 (I, J) :X (I, J, K) ) =D (K) ) ; 

©FOR (TRANSPORT (I, J,K) :X(I, J,K) >=0) ; 

©FOR (WAREHOUSE (J) : @SUM (DUMMY3 ( I , K) :X (I, J, K) ) <=CAP ( J) *Y ( J) ) 
@FOR (WAREHOUSE (J) :@SUM(DUMMY3 (I,K) :X(I, J,K) ) <=Y ( J) ) ; 

@FOR (WAREHOUSE ( J) :@BIN (Y ( J) ) ) ; 


FORMULATION S1V1.1 (B) 


SETS : 

SUPPLY/ 1. . 20/ : S ; 

WAREHOUSE/1. . 20/ : CAP, F, Y; 

DEMAND/1. .20/ :D; 

TRANSPORT ( SUPPLY , DEMAND , WAREHOUSE ) :C,X; 
DUMMY1 (SUPPLY, DEMAND) ; 

DUMMY 2 (WAREHOUSE , DEMAND ) ; 

DUMMY 3 (SUPPLY, WAREHOUSE) ; 

ENDSETS 


DATA: 


ENDDATA 


MIN=@SUM (TRANSPORT (I , J, K) :C (I, J, K) *X (I, J, K) ) + 
@SUM (WAREHOUSE (J) :F(J) *Y (J) ) ; 



! SUBJECT TO 


@SUM (TRANSPORT (I , J,K) :X(I, J,K) )=1; 

@FOR (SUPPLY (I) : @SUM (DUMMY2 ( J, K) :X (I, J,K) ) <=S (I) ) ; 

@FOR (DEMAND (K) :@SUM(DUMMY3 (I, J) :X (I , J, K) ) =D (K) ) ; 

@FOR (TRANSPORT (I, J,K) :X(I, J,K) >=0) ; 

@FOR (WAREHOUSE (J) :@SUM (DUMMY3 (I , K) :X (I , J, K) ) <=CAP ( J) *Y ( J) ) ; 
@FOR (WAREHOUSE (J) :@SUM(DUMMY3 (I , K) :X (I, J, K) )<=Y(J) } ; 

©FOR (WAREHOUSE (J) : Y ( J) >=0) ; 


FORMULATION S1V1.2(B) 


SETS : 

SUPPLY/1. .20/ : S ; 

WAREHOUSE/1 . . 20/ : CAP, F, Y; 

DEMAND/1. .20/ : D ; 

TRANS PORT ( SUPPLY , DEMAND , WAREHOUSE ) : C , X ; 
DUMMY1 (SUPPLY, DEMAND) ; 

DUMMY2 (WAREHOUSE, DEMAND) ; 

DUMMY3 (SUPPLY, WAREHOUSE) ; 

ENDSETS 


DATA: 


M=100000 ; 

ENDDATA 

MIN=@SUM (TRANSPORT (I, J, K) :C(I, J, K) *X(I, J,K) ) + 

©SUM (WAREHOUSE ( J) : F ( J) *Y ( J) ) ; 

'SUBJECT TO 

@SUM (TRANSPORT (I, J, K) :X(I, J, K) ) =1; 

@FOR (SUPPLY (I) : @SUM (DUMMY2 ( J, K) :X(I, J,K) ) <=S (I) ) ; 

©FOR (DEMAND (K) :@SUM(DUMMY3 (I, J) :X(I, J,K) ) =D (K) ) ; 

@FOR (TRANSPORT ( I, J,K) :X(I, J,K) >=0) ; 

@FOR (WAREHOUSE (J) :@SUM(DUMMY1 (I,K) :X(I, J, K) ) <=CAP ( J) *Y (J) ) ; 

@FOR (WAREHOUSE (J) :@SUM(DUMMY1 (I,K) :X(I, J,K) ) >=M*Y (J) -M) ; 
@FOR (WAREHOUSE (J) :@SUM (DUMMY1 (I , K) :X(I, J,K) ) >=-M*Y(J) ) ; 
@FOR (WAREHOUSE (J) :@SUM (DUMMY1 (I, K) :X(I, J,K) ) <=M*Y (J) ) ; 

@FOR (WAREHOUSE (J) :Y(J)>=0) ; 


FORMULATION S1V1.3(B) 

SETS : 

SUPPLY/1. .20/ : S ; 

WAREHOUSE/1. .20/ :CAP,F,Y; 

DEMAND/1. .20/ :D; 

TRANSPORT (SUPPLY , DEMAND , WAREHOUSE) :C,X; 
DUMMY1 (SUPPLY, DEMAND) ; 

DUMMY2 (WAREHOUSE, DEMAND) ; 



DUMMY 3 (SUPPLY, WAREHOUSE) ; 
ENDSETS 


DATA: 


ENDDATA 

MIN=@SUM (TRANSPORT (I ,J,K) : C (I, J, K) *X (I, J, K) ) + 

@SUM (WAREHOUSE ( J) : F ( J) *Y ( J) ) ; 

! SUBJECT TO 

@SUM (TRANSPORT (I, J,K) :X(I, J,K) ) =1; 

@FOR (SUPPLY (I) : @SUM (DUMMY2 (J,K) :X (I, J, K) ) <=S (I) ) ; 

@FOR (DEMAND (K) :@SUM(DUMMY3 (I, J) :X (I, J,K) ) =D (K) ) ; 

@FOR (TRANSPORT ( I, J,K) :X(I, J, K) >=0) ; 

@FOR (WAREHOUSE ( J) :@SUM ( DUMMY 1 (I, K) :X (I, J, K) ) <=CAP (J) *Y (J) ) ; 

@FOR (DEMAND (K) : @FOR (WAREHOUSE (J) :@SUM (SUPPLY (I) :X (I , J, K) ) <=D (K) *Y(J))) 

@FOR (WAREHOUSE (J).:Y(J) >=0) ; 


FORMULATION S1V1.4(B) 


SETS : 

SUPPLY/1. .20/ :S; 

WAREHOUSE/1. . 20/ : CAP , F, Y; 

DEMAND/1. .20/ :D; 

TRANSPORT (SUPPLY, DEMAND, WAREHOUSE) : C,X; 
DUMMY1 (SUPPLY, DEMAND) ; 

DUMMY2 (WAREHOUSE, DEMAND) ; 

DUMMY3 (SUPPLY, WAREHOUSE) ; 

ENDSETS 


DATA: 


M=100000 ; 

ENDDATA 

MIN=@SUM (TRANSPORT ( I, J,K) :C(I,J,K) *X(I,J,K) ) + 

@SUM (WAREHOUSE ( J) : F ( J) *Y ( J) ) ; 

! SUBJECT TO 

@SUM (TRANS PORT ( I, J,K) :X(I, J,K) ) =1; 

@FOR (SUPPLY (I) : @SUM (DUMMY 2 (J,K) :X (I, J,K) ) < = S (I) ) ; 

@FOR (DEMAND (K) :@SUM(DUMMY3 (I, J) :X(I, J,K) ) =D (K) ) ; 

@FOR (TRANSPORT (I, J,K) :X(I, J,K) >=0) ; 

@FOR (WAREHOUSE (J) : @SUM (DUMMY1 ( I , K) :X(I, J,K) ) <=CAP ( J) *Y ( J) ) ; 

@FOR (DEMAND (K) :@FOR (WAREHOUSE (J) :@SUM (SUPPLY(I) :X (I, J, K) ) < = - 
M*Y ( J) +M+D (K) ) ) ; 

@FOR (DEMAND (K) :@FOR (WAREHOUSE (J) :@SUM (SUPPLY (I) :X (I, J, K) ) > = -M*Y (J) ) ) ; 
@FOR (DEMAND (K) :@FOR (WAREHOUSE ( J) :@SUM (SUPPLY (I) :X (I, J, K) ) <=M*Y ( J) ) ) ; 
@FOR (WAREHOUSE ( J) :Y (J) >=0) ; 



FORMULATION S1V1.5(B) 


SETS: 

SUPPLY/ 1. . 20/ : S ; 

WAREHOUSE/1. . 2 0/ : CAP, F, Y; 

DEMAND/1. .20/ :D; 

TRANSPORT ( SUPPLY , DEMAND , WAREHOUSE ) : C , X ; 
DUMMY 1 (SUPPLY, DEMAND) ; 

DUMMY2 (WAREHOUSE, DEMAND) ; 

DUMMY3 (SUPPLY, WAREHOUSE) ; 

ENDSETS 


DATA: 


ENDDATA 

MIN=@SUM (TRANSPORT (I, J,K) : C (I , J, K) *X (I, J, K) ) + 
@SUM (WAREHOUSE ( J) : F ( J) *Y ( J) ) ; 


! SUBJECT TO 

©SUM (TRANSPORT ( I, J,K) :X(I, J,K) )=1; 

@FOR (SUPPLY (I) :@SUM (DUMMY2 (J,K) :X(I, J,K) ) <=S (I) ) ; 

@FOR (DEMAND (K) :@SUM(DUMMY3 (I, J) :X (I , J, K) ) =D (K) ) ; 

@FOR (TRANSPORT (I, J, K) :X(I, J, K) >=0) ; 

@FOR (WAREHOUSE (J) :@SUM(DUMMY1 (I,K) :X(I, J, K) ) <=CAP (J) *Y (J) ) ; 

@FOR (DEMAND (K) :@FOR (WAREHOUSE (J) :@SUM (SUPPLY (I) :X (I, J, K) ) <=D (K) *Y ( J) ) ) 
@FOR (WAREHOUSE (J) :@FOR (SUPPLY (I ) : @SUM (DEMAND (K) :X (I, J,K) ) <=S (I) *Y(J) ) ) 

@FOR (WAREHOUSE ( J) : Y ( J) >=0) ; 


FORMULATION S1V1.6 (B) 


SETS: 

SUPPLY/1. .20/ : S ; 

WAREHOUSE/1. .20/ :CAP,F,Y; 

DEMAND/1. .20/ :D; 

TRANSPORT (SUPPLY, DEMAND, WAREHOUSE) :C,X; 
DUMMY1 (SUPPLY, DEMAND) ; 

DUMMY2 (WAREHOUSE, DEMAND) ; 

DUMMY 3 (SUPPLY, WAREHOUSE) ; 

ENDSETS 


DATA: 

M= 100000; 

ENDDATA 

MIN=@SUM (TRANSPORT (I, J, K) : C (I , J, K) *X (I , J, K) ) + 
@SUM (WAREHOUSE (J) :F (J) *Y (J) ) ; 


! SUBJECT TO 



©SUM (TRANSPORT (I,J,K) :X(I,J,K))=1? 

@FOR (SUPPLY (I) : @SUM (DUMMY2 (J,K) :X(I, J,K) ) <=S(I) ) ; 

@FOR (DEMAND (K) :@SUM(DUMMY3 (I, J) :X(I,J,K) )=D(K) ) ; 

@FOR (TRANSPORT (I ,J,K) :X(I, J,K) >=0) ; 

@FOR (WAREHOUSE (J) : @SUM (DUMMY1 ( I , K) :X(I,J,K))<=CAP(J)*Y(J)); 

@FOR (DEMAND (K) : @FOR (WAREHOUSE (J) : @SUM (SUPPLY (I) :X (I , J, K) ) < = - 
M*Y ( J) +M+D (K) ) ) ; 

@FOR (DEMAND (K) :@FOR (WAREHOUSE (J) :@SUM (SUPPLY (I) :X (I , J, K) ) >=-M*Y ( J) ) ) ; 
@FOR (DEMAND (K) : @FOR (WAREHOUSE (J) : @SUM ( SUPPLY ( I ) :X (I , J, K) ) <=M*Y ( J) ) ) ; 

@FOR (WAREHOUSE (J) :@FOR (SUPPLY (I) :@SUM (DEMAND (K) :X (I , J, K) ) <=- 
M*Y ( J) +M+S (I) ) ) ; 

©FOR (WAREHOUSE (J) :@FOR (SUPPLY (I) : @SUM (DEMAND (K) :X (I , J, K) ) >=-M*Y ( J) ) ) ,- 
©FOR (WAREHOUSE (J) :@FOR (SUPPLY (I) : ©SUM (DEMAND (K) :X(I,J,K) )<=M*Y(J) )) ; 

@FOR (WAREHOUSE (J) :Y(J) >=0) ; 


FORMULATION S1V1.7 (B) 

SETS: 

SUPPLY/1. . 20/ : S ; 

WAREHOUSE/1. . 20/ : CAP , F, Y; 

DEMAND/1. .20/ :D; 

TRANSPORT (SUPPLY, DEMAND, WAREHOUSE) : C, X; 
DUMMY 1 (SUPPLY, DEMAND) ; 

DUMMY2 (WAREHOUSE, DEMAND) ; 

DUMMY3 (SUPPLY, WAREHOUSE) ; 

ENDSETS 


DATA: 


ENDDATA 


MIN=@SUM ( TRANSPORT ( I, J,K) :C (I , J, K) *X (I , J, K) ) + 
@SUM (WAREHOUSE ( J) : F ( J) *Y ( J) ) ; 


! SUBJECT TO 

©SUM (TRANSPORT (I, J,K) :X(I, J,K) )=1; 

@FOR (SUPPLY (I) :@SUM (DUMMY 2 (J, K) :X (I, J, K) ) <=S (I) ) ; 

©FOR (DEMAND (K) :@SUM (DUMMY 3 (I, J) :X (I , J, K) ) =D (K) ) ? 

©FOR (TRANSPORT (I, J,K) :X(I, J,K) >=0) ; 

@FOR (WAREHOUSE (J) :@SUM(DUMMY3 (I, K) :X (I, J,K) ) <=CAP (J) *Y (J) ) ; 
@FOR (WAREHOUSE ( J) : Y ( J) >=0) ; 



STYLE 1 VARIATION 2 


Style 1 Variation 2 has formulations same as the formulations for Style 1 Variation 1 
except that the constraint (in LINGO) 

@FOR(WAREHOUSE(J):@SUM(DUMMY3(I,K):X(I,J,K))<=CAP(J)*Y(J)); 

Is replaced by constraint 

@FOR(WAREHOUSE(J):@SUM(DUMMY3(I,K):X(I,J,K))<=CAP(J)); 

Rest of the structure is identically same. 


STYLE 2 VARIATION 1 


FORMULATION S2Vl. 1 (A) 


SETS : 

SUPPLY/1. .20/ :S; 

WAREHOUSE/1. . 20/ :CAP, F, Y; 
DEMAND/1. .20/:D; 

SUPP_WAR ( SUPPLY , WAREHOUSE ) : Cl , XI ; 
WAR_DEM (WAREHOUSE , DEMAND ) : C2 , X2 ; 

ENDSETS 


DATA: 


ENDDATA 

MIN=@SUM(SUPP_WAR(I, J) :C1(I,J) *X1(I,J) ) + 

@SUM(WARJ0EM(J,K) :C2 (J,K) *X2 (J,K) ) + 

@SUM (WAREHOUSE (J) :F(J) *Y(J) ) ; 

! SUBJECT TO 

@FOR (SUPPLY (I) :@SUM (WAREHOUSE (J) :X1 (I, J) ) <=S (I) ) ; 

@FOR (DEMAND (K) :@SUM (WAREHOUSE ( J) :X2 ( J, K) ) >=D (K) ) ; 

©SUM (WAREHOUSE ( I, J) :X1(I,J))=1; 

@SUM(WAR_DEM(J,K) :X2 (J,K))=1; 

©FOR ('WAREHOUSE (J) :@FOR (SUPPLY (I) :X1 (I, J) >=0) ) ; 

@FOR (DEMAND (K) : @FOR (WAREHOUSE (J) :X2 (J,K)>=0) ) ; 

@FOR (WAREHOUSE (J) : @SUM ( SUPPLY ( I ) :X1 (I, J) ) =@SUM (DEMAND (K) :X2 (J,K))) 
@FOR (WAREHOUSE (J) : @SUM ( SUPPLY ( I ) :X1 (I, J) ) <=CAP ( J) *Y ( J) ) ; 

@FOR (WAREHOUSE (J) :@SUM (SUPPLY (I) :X1 (I , J) ) <=Y { J) ) ; 

@FOR (WAREHOUSE ( J) :@BIN (Y ( J) ) ) ; 


FORMULATION S2V1.1(B) 

SETS: " 

SUPPLY/1. .20/ :S; 

WAREHOUSE/1. . 2 0/ : CAP , F, Y; 
DEMAND/1. .20/ :D; 

SUPP_WAR (SUPPLY, WAREHOUSE) : Cl , XI ; 
WAR_DEM (WAREHOUSE, DEMAND) :C2,X2; 

ENDSETS 

DATA: 


ENDDATA 

MIN=@SUM(SUPP_WAR(I, J) :C1 (I, J) *X1 (I, J) ) + 
@SUM(WAR_DEM(J,K) :C2 (J,K) *X2 (J,K) ) + 
@SUM (WAREHOUSE ( J) : F ( J) *Y ( J) ) ; 


! SUBJECT TO 



@FOR (SUPPLY (I) :@SUM (WAREHOUSE (J) :X1 (I, J} ) <=S (I) ) ; 

@FOR (DEMAND (K) : @SUM (WAREHOUSE (J) :X2 (J,K))=D(K)); 

@SUM(SUPP_WAR(I, J) :X1 (I , J) ) =1 ; 

@SUM (WAR_DEM ( J, K) :X2 (J,K))=1; 

©FOR (WAREHOUSE (J) : @FOR ( SUPPLY ( I ) :X1 (I, J) >=0) ) ; 

@FOR (DEMAND (K) :@FOR (WAREHOUSE (J) :X2 (J, K)>=0)); 

@FOR (WAREHOUSE (J) : @SUM (SUPPLY (I) :X1 (I , J) ) =@SUM (DEMAND (K) :X2(J,K))) 
@FOR (WAREHOUSE (J) : @SUM (SUPPLY (I) :X1 (I , J) ) <=CAP ( J) *Y ( J) ) ; 

@FOR (WAREHOUSE (J) :@SUM (SUPPLY (I) :X1 (I, J) ) <=Y ( J) ) ; 

@FOR (WAREHOUSE (J) :Y(J) >=0) ; 


FORMULATION S2V1.2(B) 

SETS: 

SUPPLY/1. .20/ :S; 

WAREHOUSE/1. . 20/ : CAP , F , Y; 
DEMAND/1. . 2 0 / : D ; 

SUPP_WAR ( SUPPLY , WAREHOUSE ) : Cl , XI ; 
WAR_DEM ( WAREHOUSE , DEMAND ) : C2 , X2 ; 

ENDSETS 

DATA: 


M=100000; 

ENDDATA 

MIN=@SUM (SUPP_WAR ( I , J) : Cl ( I , J) *X1 ( I , J) ) + 

@SUM(WAR_DEM(J,K) :C2 ( J, K) *X2 (J,K) ) + 

@SUM (WAREHOUSE (J) :F (J) *Y (J) ) ; 

! SUBJECT TO 

@FOR (SUPPLY (I) :@SUM (WAREHOUSE (J) :X1 (I, J) ) <=S (I) ) ; 

@FOR (DEMAND (K) : @SUM (WAREHOUSE (J) :X2(J,K))=D(K)); 

@SUM(SUPP_WAR(I, J) :X1 (I, J) ) =1; 

@SUM (WAR_DEM ( J, K) :X2(J,K))=1; 

@FOR (WAREHOUSE (J) : @FOR (SUPPLY ( I) :X1 (I, J) >=0) ) ; 

@FOR (DEMAND (K) : @FOR (WAREHOUSE (J) :X2(J,K)>=0)); 

@ FOR (WAREHOUSE (J) : @SUM (SUPPLY ( I ) :X1 (I, J) ) =@SUM (DEMAND (K) :X2 (J,K) ) ) 
@FOR (WAREHOUSE (J) :@SUM (SUPPLY (I) :X1 (I, J) ) <=CAP ( J) *Y ( J) ) ; 

@FOR (WAREHOUSE (J) : @SUM (SUPPLY ( I ) :X1(I,J))>=M*Y(J)-M); 

@FOR (WAREHOUSE (J) : ©SUM (SUPPLY ( I ) : XI (I , J) )> = -M*Y(J) } ; 

@FOR (WAREHOUSE (J) :@SUM (SUPPLY (I) :X1 (I , J) )<=M*Y(J) ) ; 

@FOR (WAREHOUSE (J) :Y(J) >=0) ; 


FORMULATION S2V1.3(B) 


SETS : 

SUPPLY/ 1. .20/ :S; 



WAREHOUSE/1. . 20/ : CAP, F, Y; 
DEMAND/1. .20/ : D ; 

SUPP_WAR (SUPPLY, WAREHOUSE) :C1,X1; 
WAR_DEM (WAREHOUSE , DEMAND) : C2 , X2 ; 

ENDSETS 

DATA: 

ENDDATA 


MIN=@SUM ( SUPP_WAR ( I , J) : Cl ( I , J) *X1 (I , J) ) + 
@SUM(WAR_DEM(J,K) :C2 ( J, K) *X2 (J,K) ) + 
@SUM (WAREHOUSE ( J) : F ( J) *Y ( J) ) ; 


! SUBJECT TO 

©FOR (SUPPLY (I) :@SUM (WAREHOUSE (J) :X1 (I, J) ) <=S (I) ) ; 

@FOR (DEMAND (K) : @SUM (WAREHOUSE (J) :X2 ( J, K) ) =D (K) ) ; 

@SUM ( SUPP_WAR ( I , J) :X1 (I, J) )=1; 

@SUM(WAR_DEM(J,K) :X2(J,K))=1; 

@FOR (WAREHOUSE (J) :@FOR (SUPPLY (I) :X1 (I, J) >=0) ) ; 

@FOR (DEMAND (K) : @FOR (WAREHOUSE (J) :X2(J,K) >=0)); 

@FOR (WAREHOUSE (J) : @SUM (SUPPLY (I) :X1 (I , J) ) =@SUM (DEMAND (K) :X2 (J,K))) 
@FOR (WAREHOUSE (J) :@SUM (SUPPLY (I) :X1 (I, J) ) <=CAP (J) *Y (J) ) ; 
@FOR(WAR_DEM(J,K) :X2 ( J, K) <=Y ( J) *D (K) ) ; 

@FOR (WAREHOUSE (J) :Y(J) >=0) ; 


FORMULATION S2V1.4(B) 


SETS: 

SUPPLY/1. .20/ :S; 

WAREHOUSE/1. . 20/ : CAP, F , Y; 

DEMAND/1. .20/ :D; 

SUPP_WAR ( SUPPLY , WAREHOUSE ) : Cl , XI ; 

WAR_DEM (WAREHOUSE , DEMAND) : C2 , X2 ; 

ENDSETS 

DATA: 

M=100000 ; 

ENDDATA 

MIN=@SUM (SUPP_WAR ( I , J) :C1(I, J) *X1 (I, J) ) + 

@SUM (WAR_DEM ( J, K) :C2 (J,K) *X2 (J,K) ) + 

@SUM (WAREHOUSE ( J) : F ( J) *Y ( J) ) ; 

! SUBJECT TO 

@FOR (SUPPLY (I) :@SUM (WAREHOUSE (J) :X1 (I, J) ) <=S (I) ) ; 


@FOR (DEMAND (K) :@SUM (WAREHOUSE ( J) :X2 ( J, K) ) =D (K) ) ; 

@SUM(SUPP_WAR(I, J) :X1 (I, J) ) =1; 

@SUM(WAR_DEM(J,K) :X2 (J,K) )=1; 

@FOR (WAREHOUSE (J) : ©FOR (SUPPLY (I ) :X1 (I, J) >=0) ) ; 

@FOR (DEMAND (K) :@FOR (WAREHOUSE ( J) :X2 (J, K) >=0) ) ; 

©FOR (WAREHOUSE (J) :@SUM (SUPPLY (I) :X1 (I , J) ) =@SUM (DEMAND (K) :X2 (J,K))) 
@FOR (WAREHOUSE (J) :@SUM (SUPPLY (I) :X1 (I , J) ) <=CAP ( J) *Y ( J) ) ; 
@FOR(WAR_DEM(J,K) :X2 ( J, K) < = -M*Y ( J) +M+D (K) ) ; 

@FOR (WAR_DEM ( J, K) :X2 ( J, K) <=M*Y ( J) ) ; 

@FOR(WAR_DEM(J,K) :X2 (J,K) >=-M*Y(J) ) ; 

@FOR (WAREHOUSE ( J) : Y ( J) >=0) ; 


FORMULATIONS S2V1. (B) 


SETS: 

SUPPLY/1. .20/ :S; 

WAREHOUSE/1. .20/: CAP ,F,Y; 
DEMAND/1. .2 0/-.D; 

SUPP_WAR( SUPPLY, WAREHOUSE) :C1,X1; 
WAR_DEM (WAREHOUSE , DEMAND) : C2 , X2 ; 

ENDSETS 


DATA: 


ENDDATA 


MIN=@SUM (SUPP_WAR ( I , J) : Cl ( I , J) *X1 ( I , J) ) + 

@SUM(WAR_DEM(J,K) :C2 (J,K) *X2 (J,K) ) + 

@SUM (WAREHOUSE ( J) : F ( J) *Y ( J) ) ; 

! SUBJECT TO 

@FOR (SUPPLY (I ) :@SUM (WAREHOUSE (J) :X1 (I, J) )<=S(I)) ; 

@FOR (DEMAND (K) : @SUM (WAREHOUSE (J) :X2 (J, K) ) =D (K) ) ; 

@SUM(SUPP_WAR(I, J) : XI ( I , J) ) =1 ; 

@SUM (WAR_DEM ( J, K) :X2(J,K))=1; 

©FOR (WAREHOUSE (J) :@FOR (SUPPLY (I) :X1 (I, J) >=0) ) ; 

@FOR (DEMAND (K) : @FOR (WAREHOUSE (J) :X2(J,K)> = 0) ) ; 

@FOR (WAREHOUSE (J) :@SUM (SUPPLY (I) :X1 (I , J) ) =@SUM (DEMAND (K) :X2(J,K) ) ) 
©FOR (WAREHOUSE (J) : @SUM (SUPPLY ( I ) :X1 (I, J) ) <=CAP ( J) *Y ( J) ) ; 

@FOR (WAR_DEM ( J, K) :X2 ( J, K) <=Y ( J) *D (K) ) ; 

@FOR(SUPP_WAR(I, J) : XI ( I , J) <=Y ( J) *S(I) ) ; 

@FOR (WAREHOUSE (J) :Y(J)>=0) ; 


FORMULATION S2V1 . 6 (B) 


SETS : 

SUPPLY/1. .20/ :S; 
WAREHOUSE/1. .20/:CAP,F,Y; 
DEMAND/1. .20/ :D; 



SUPP_WAR ( SUPPLY , WAREHOUSE ) : Cl , XI ; 
WAR_DEM ( WAREHOUS E , DEMAND ) : C2 , X2 ; 

ENDSETS 


DATA: 


M=100000 ; 

ENDDATA 

MIN=@SUM ( SUPP_WAR ( I , J) : Cl ( I , J) *X1 ( I , J) ) + 

@SUM (WAR_DEM ( J, K) : C2 ( J, K) *X2 ( J, K) ) + 

@SUM (WAREHOUSE ( J) : F ( J) *Y ( J) ) ; 

! SUBJECT TO 

@FOR (SUPPLY (I) :@SUM (WAREHOUSE (J) :X1 (I, J) ) <=S (I) ) ; 

@FOR (DEMAND (K) :@SUM (WAREHOUSE ( J) :X2 ( J, K) ) =D (K) ) ; 

@SUM (SUPP_WAR ( I , J) :X1 (I , J) ) =1 ; 

@SUM(WAR_DEM(J,K) :X2 (J,K))=1; 

@FOR (WAREHOUSE (J) :@FOR (SUPPLY (I) :X1 (I, J) >=0) ) ; 

@FOR (DEMAND (K) :@FOR (WAREHOUSE ( J) :X2 ( J, K) >=0) ) ; 

@FOR (WAREHOUSE (J) : @SUM (SUPPLY (1) :X1 (I, J) ) =@SUM (DEMAND (K) :X2(J,K))); 
@FOR (WAREHOUSE (J) :@SUM (SUPPLY (I) :X1 (I, J) ) <=CAP ( J) *Y (J) ) ; 

@FOR (WAR_DEM ( J, K) :X2 ( J, K) <=-M*Y ( J) +M+D (K) ) ; 

@FOR (WAR_DEM ( J, K) :X2 ( J, K) <=M*Y ( J) ) ; 

@FOR (WAR_DEM ( J, K) :X2 ( J, K) >=-M*Y ( J) ) ; 

@FOR(SUPP_WAR(I, J) :X1(I, J) <=-M*Y(J) +M+S (I) ) ; 

@FOR(SUPP_WAR(I, J) :X1(I, J)<=M*Y(J) ) ; 

@FOR ( SUPP_WAR ( I , J) :X2 (I, J) >=-M*Y(J) ) ; 


@FOR (WAREHOUSE (J) :Y(J) >=0) ; 


FORMULATION S2V1.7(B) 


SETS : 

SUPPLY/1. .20/ :S; 

WAREHOUSE/1. . 20/ : CAP , F , Y; 
DEMAND/1. .20/ :D; 

SUPP_WAR ( SUPPLY , WAREHOUSE ) : Cl , XI ; 
WAR_DEM (WAREHOUSE , DEMAND) : C2 , X-2 ; 

ENDSETS 


DATA: 


ENDDATA 


MIN=@SUM (SUPP WAR ( I , J) :C1 (I, J) *X1 (I, J) ) + 



@SUM(WAR_DEM(J,K) :C2 ( J, K) *X2 (J,K) ) + 
@SUM (WAREHOUSE ( J) : F ( J) *Y ( J) ) ; 


! SUBJECT TO 

@FOR (SUPPLY (I) :@SUM (WAREHOUSE (J) :X1 (I, J) ) <=S (I) ) ; 

@FOR (DEMAND (K) :@SUM (WAREHOUSE (J) :X2 (J, K) ) =D (K) ); 

@SUM(SUPP_WAR(I, J) :X1 ( I , J) ) =1 ; 

@SUM(WAR_DEM(J,K) :X2 (J,K))=1; 

@FOR (WAREHOUSE (J) : @FOR (SUPPLY ( I) :X1(I,J) >=0)) ? 

@FOR (DEMAND (K) :@FOR (WAREHOUSE ( J) :X2 ( J, K) >=0) ) ; 

@FOR (WAREHOUSE (J) :@SUM (SUPPLY (I) :X1 (I, J) ) =@SUM (DEMAND (K) :X2(J,K) ) ) ; 
@FOR (WAREHOUSE (J) :@SUM (SUPPLY (I) :X1 (I, J) ) <=CAP ( J) *Y ( J) ) ; 

@FOR (WAREHOUSE ( J) : Y ( J) >=0 ) ; 


STYLE 2 VARIATION 2 


Style 1 Variation 2 has formulations same as the formulations for Style 1 Variation 1 
except that the constraint (in LINGO) 

@FOR (WAREHOUSE (J) : ©SUM (SUPPLY (I) :X1 (I , J) ) <=CAP ( J) *Y ( J) ) ; 

Is replaced by constraint 

@FOR (WAREHOUSE (J) :@SUM (SUPPLY (I) :X1 (I , J) ) <=CAP ( J) ) ; 

Rest of the structure is identically same. 



